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BIOLOGICAL  EFFECTS  OF  MILLIMETER-WAVE  IRRADIATION 


INTRODUCTION 

The  increased  use  of  millimeter  waves  for  civilian  and  military  purposes 
has  generated  great  interest  in  the  possible  health  hazards  of  millimeter-wave 
radiation  at  thermal  and  nonthermal  levels  [1].  Also,  the  presence  of  several 
reports  indicating  that  biological  effects  could  be"  induced  by  low-level 
radiation  (<5  mW/cm  )  has  stimulated  interest  in  the  more  basic  question  of 
the  mechanisms  of  interaction  of  millimeter  waves  with  biological  systems. 

This  project  has  approached  the  question  of  the  biological  effects  of 
millimeter  waves  on  several  fronts:  independent  verification  of  the  reported 
frequency-sensitive  growth  rates  of  yeast  cultures  to  millimeter-wave 
radiation  [2-4],  Raman  spectroscopic  studies  of  bacterial  cell  suspensions  and 
of  a  model  of  cell  membranes  [liposomes],  determination  of  dielectric 
.properties  of  biologically  relevant  samples,  and  experiments  to  determine  the 
presence  of  surface  acoustic  waves  in  polystyrene  microspheres  and  liposomes. 

EFFECTS  OF  MILLIMETER-WAVE  IRRADIATION  ON 
GROWTH  OF  SACCHAROMYCES  CEREVISIAE 

The  growth  of  Saccharomyces  cerevisiae  cultures  has  been  reported  [2-4] 
to  be  strongly  influenced  by  millimeter-wave  irradiation  around  41.7  GHz.  In 
particular,  a  change  of  only  8  MHz  in  the  irradiation  frequency  was  reported 
to  be  sufficient  to  alter  significantly  the  magnitude  of  the  effect. 

The  need  to  verify  such  extreme  sensitivity  of  yeast  cells  has  motivated 
this  section  of  the  project.  A  new  irradiation  chamber  was  designed  and  built 
to  allow  simultaneous  irradiation  and  sham  irradiation  of  recirculated  cell 
suspensions.  The  temperature  difference  was  <  0.01°C  and  the  growth  was 
measured  optically.  A  klystron,  phase-locked  by  an  EIP-578  microwave  counter, 
was  used  to  generate  millimeter  waves  stable  to  within  ±  50  Hz  for  the  4  h  of 
irradiation  and  the  frequency  was  measured  with  up  to  11  significant  digits. 
At  least  three  experiments  were  done  at  each  frequency.  The  irradiation 
frequencies  were  in  the  ranges,  between  41.650  and  41.798  GHz,  in  which  the 


largest  -effects  were  previously  ^eoortec!  [2-1*].  Tr,  these  regions  the  frequen¬ 
cies  were  2  MHz  apart. 

No  difference  larger  than  ±  U%  was  detected  in  the  yeast  growth  rate  at 
any  of  the  irradiation  frequencies.  Furthermore,  these  effects  were  not  sig¬ 
nificant  at  the  95'-  confidence  limit.  Plate  counts  performed  at  the  end  of 
each  experiment  were  in  fair  agreement  with  the  growth  data.  A  complete 
report  on  these  experiments  has  been  published  [5]  and  its  reprint  is  attached 
as  Appendix  A. 

RAMAN  SPECTROSCOPY  OF  LIPOSOMES 
EXPOSED  TO  MILLIMETER  WAVES 

Cultures  of  Bacillus  negater  turn  were  reported  to  show  Raman  spect-a  de¬ 
pendent  on  the  stage  of  their  life  cycle  [6,  7].  A  preliminary  study  was  con¬ 
ducted  to  investigate  the  speqtral  features  of  Raman  spectra  of  megater ium 
cells  after  heat  and  cold  shocks,  as  indicated  in  references  5  and  7.  No 
biologically  related  features  were  detected  in  the  Raman  spectra  at  any  stage 
of  the  cell  cycle.  These  results  have  been  published  [8,  9],  and  their  re¬ 
paints  are  attached  as  Appendixes  5  and  C. 

Eve-1  if  Raman  spectra  had  be»r.  detected  for  cultures  of  bacterial  cells, 
i f  would  hav°  been  rather  difficult  to  ascribe  their  origin  to  a  particular 
structure  of  molecular  process.  We,  therefore,  decided  to  use  another  ap¬ 
proach  to  study  the  possible  mechanism  of  interaction  of  millimeter  waves  with 
biological  systems.  A  relatively  simple  and  well-characterized  model  of 
cytoplasmic  membranes  has  been  used;  and  Raman  spectroscopy  was  used  to  prove, 
concurrently  with  the  irradiation,  the  ability  of  millimeter  waves  to  induce 
conformational  changes  in  lipid  bilayers.  Experiments  were  performed  on 
d ipalmitoy l-phosphat idy lchol ine  (DPPC)  liposomes  below  and  above  the  transi¬ 
tion  temperature,  ^i°C,  of  DPPC.  As  expected,  a  temoerature  increase  induced 
spectral  changes,  while  no  differences  were  detected  upon  millimeter  wave 
irradiation,  a  comolete  r^Dcrt  on  these  experiments  has  b«en  published  in  ‘“n 
Journal  of  Applied  Physics  I i 0 ] .  A  reprint  coDy  of  reference  -0  is  attached 
as  Amend ix  D. 
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BFILLOUIN  SPECTROSCOPY  OF  LIPOSOMES 
AND  MICROSPHERES 


It  is  possible  that  surface  acoustic  waves  may  be  supported  on 
tnicrospheres  and  also  on  liposomes.  To  test  this  hypothesis,  preliminary 
experiments  have  been  performed  with  Brillouin  spectroscopy  using  a  Fabry- 
Perot  interferometer  in  single  and  five-pass  configurations.  The  spectra  of 
CH C.£i*  and  of  water  are  shown  in  Figures  1  and  2.  These  spectra  were  obtained 
by  accumulating  100  spectra  on  1024  channels  with  2-ms  accumulation/channe 1 . 
The  laser  source  was  set  at  488.8  nm  and  an  etalon  was  placed  in  the  laser. 
Experiments  are  presently  underway  to  determine  the  presence  of  surface 
acoustic  waves  in  microspheres  and  liposomes. 

EXPERIMENTAL  DETERMINATION  OF  COMPLEX 
PERMITTIVITIES  OF  BIOLOGICAL  SAMPLES 

We  have  developed  a  computer-controlled  measurement  system  operating  at 
frequencies  from  26.5  to  60  GHz.  The  measurement  system  is  schematically 
illustrated  in  Figure  3.  The  key  components  of  the  measurement  system  are  a 
precision  slotted  line,  a  diode  detector  for  standing-wave  measurements,  a 
frequency  counter  EIP-578,  and  HP  86B  computer  for  overall  automation.  The 
traveling  probe  of  the  slotted  line  is  moved  in  steps  of  0.3  mn  by  a  computer- 
controlled  stepping  motor.  A  large  dial  indicator  calibrated  in  0.01-mm 
increments  gives  the  position  of  the  probe.  At  each  location  of  the  slotted 
line  probe,  the  standing-wave  voltage  is  measured  by  the  diode  detector 
(sensitivity:  -  60  dBm)  via  SWR  meter,  A/D  converter  and  HP  86B  computer.  By 
averaging  10  or  more  repetitive  measurements  of  standing-wave  voltages,  the 
extreme  variation  was  limited  to  ±  0.06  dB.  The  standing-wave  voltages  were 
measured  at  40  equally  spaced  positions  of  the  slotted  line  probe.  The  phase 
and  magnitude  of  the  reflection  coefficient  were  obtained  by  implementing  the 
three-probe  ref lectometer  formulation  given  by  Caldecott  [11].  Between  20  to 
35  independent  values  of  the  phase  and  magnitude  of  the  reflection  coefficient 
are  obtained  due  to  different  combinations  of  the  standing-wave  voltages.  The 
accuracy  of  the  reflection  coefficient  measurements  is  improved  significantly 
by  averaging  these  reflection  coefficient  values. 

Using  the  measurement  system  of  Figure  3,  the  complex  permittivities  were 
measured  for  several  samples;  e.g.,  10%  saline  solution,  human  blood,  20% 


bovine  serum  albumin  (BSA)  solution  and,  preparations  of  Saccharomyces 
cerevisiae  with  densities  of  10  and  10  cells/m£.  The  measurements  were  not 
performed  in  the  frequency  range  of  35  to  40  GHz  because  of  very  low  output 
from  Che  tunable  probe  of  Che  slotted  line.  Also,  in  Che  frequency  range  of 
46  to  50  GHz,  the  complex  permittivity  data  were  not  obtained  because  of 
problems  with  the  stability  of  IMPATT  diode  oscillator. 

A  modified  "infinite"  sample  method  has  been  developed  to  determine  the 

dielectric  properties  of  the  biological  samples.  In  this  method,  an 

impedance-matching  transformer  is  used  to  match  an  electrically  infinite 

biological  sample  to  air-filled  waveguide.  The  matching  transformer  consists 

of  a  X  /4  or  3A  /4  waveguide  section  of  a  suitable  low-loss  dielectric 
ge  ge 

material.  The  choice  of  the  matching  transformer  material  would  depend  on  its 
dielectric  constant,  loss  tangent,  and  machinability.  For  biological  samples, 
it  has  been  found  that  low-loss  materials  with  dielectric  constant  hetween  6 
to  10  are  suitable.  The  modified  infinite  sample  method  has  several 
advantages  over  other  methods  (e.g..  Von  Hippel's  short-circuit  technique, 
microwave  bridge  method,  etc.),  including  (1)  elimination  of  tedious  sample 
preparation  associated  with  submiliimeter  sample  lengths,  and  (2)  relatively 
easier  broadband  measurements. 

The  frequency  ranges  28  to  33  GHz,  42  to  45  GHz,  and  52  to  57  GHz  have 

been  covered  by  three  matching  transformers  of  alumina  of  lengths  X  / 4, 

g  ^ 

3A  /4  and  5A  /4,  respectively.  From  uncertainty  analysis,  it  has  been 
g£  gg 

observed  that  ±  10%  accuracy  in  complex  permittivity  of  a  biological  sample 
can  be  obtained  if  the  dielectric  constant  and  the  length  of  the  alumina 
matching  transformer  are  known  to  within  ±  1%  and  ±  5um  ,  respectively. 

Because  these  types  of  accuracies  cannot  be  realized  by  conventional  methods, 
it  is  difficult  to  characterize  the  matching  transformer  separately.  This 
problem  has  been  alleviated  by  estimating  the  matching  transformer  parameters 
from  measurements  with  liquids  of  known  dielectric  properties;  e.g.,  deionized 
water.  The  measured  complex  permittivities  of  deionized  water  having  a 
conductivity  of  200  uS/m  are  shown  in  Figure  4  with  the  best-fit  Debye  model 

(solid  line).  The  best-fit  Debye  parameters  (e  =  78.65,  e  *  4.04, and  t  = 
8.58  ps)  are  in  fairly  good  agreement  with  those  obtained  from  Schwan  et  al . 


-  6  - 
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[12].  The  estimated  dielectric  constant  of  alumina  matching  transformer  from 
deionized  water  measurements  is  9.25  which  is  in  close  agreement  with  9.15  ± 
0.16  measured  using  Von  Hippel's  shorted  line  technique  [13],  The  estimated 
lengths  of  the  matching  transformers  also  agree  closely  with  the  physically 
measured  lengths. 


The  measured  complex  permittivity  data,  along  with  the  data  from  best-fit 

Debye  model,  are  shown  in  Figures  4-9  for  10%  saline,  human  blood,  20%  BSA 

solution,  and  Saccharomyces  cerevisiae  cell  suspensions  of  densities  10^  and 
q 

10  cells/  mi.  The  Debye  equation  is  given  by 
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where,  and  are  the  static  permittivity  and  infinite-freauency  permit¬ 
tivity,  respectively.  a  and  r  represent  ionic  conductivity  and  relaxation 
time,  respectively.  The  Debye  parameters  (e  ,  e^,  t  and  a),  along  with 
standard  deviations,  are  given  in  Table  1  for  different  biological  materials 
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Figure  9.  Complex  permittivity  of  S^_  cerevisiae 
(density  ■  10  cells/mi). 


TABLE  1.  DEBYE  PARAMETERS  FOR  VARIOUS 
BIOLOGICAL  MATERIALS 


Water 

23*1 

78.7*1.1 

102  saline 
solution 

23±l 

55.9*0.9 

Human  blood 

23*1 

58.0*1.6 

Bovine  serum 
albumin  202 

C.  V. 

23*1 

58.5*1.4 

S.  cerevisiae 

24±1 

61.0*1.7 

'  109  cells/mi 

1 

S.  cerevisiae 

24*1 

78.3*2.0 

i  107  cells/mi 

i 

Huran  skin 

in  vivo 

— 

28.5*1.3 

. _ i _ i _ 

8.58*0.14 

7.10*0.13 


9.64*0.24 
7 .80±0. 19 


9.06±0.28 


9 .00±0 . 3 


10.n5±0.6 


that  have  been  tested.  The 
defined  as 


root 


mean 


square  error  (RMSE)  given  in  Table  1  is 


RMSE  - 

★  * 

where  s  .  and  e.  represent  measured  complex  permittivity  and  complex 

mi  i 

permittivity  from  the  Debye  equation  for  ith  frequency,  respectively,  and  N  is 
the  number  of  frequencies. 

MEASUREMENT  OF  THE  COMPLEX  PERMITTIVITY 
OF  HUMAN  SKIN  IN  VIVO 

The  computer-controlled  measurement  system  described  in  the  last  section 
was  used  for  computer  permittivity  measurements  of  the  human  skin  in  vivo. 
The  complex  reflection  coefficient  was  measured  by  pressure  contact  of  the 
waveguide  flange  against  an  area  of  the  body,  mostly  the  palm  of  the  hand. 
The  high  reflectivity  of  the  human  skin  was  reduced  by  using  a  quarter 

wavelength  matching  transformer  of  Teflon.  Complex  permittivities  were 
determined  from  measured  reflection  coefficients  using  a  variational  method 
iteratively  for  calculating  the  input  impedance.  A  zero-finding  technique 
algorithm  was  used  for  optimization. 

Figure  10  shows  the  real  and  imaginary  parts  of  the  permittivity  for  the 
skin  of  the  palm  along  with  the  values  obtained  from  the  best-fit  Debye 
equation.  Equation  1  for  the  Debye  relaxation  model  was  fitted  to  our  data 

for  Che  frequency  band  26.5  -  60  GHz  and  to  the  lower  frequency  data  reported 

by  Hey-Shipton  et  al.  [14],  The  Debye  parameters  for  the  skin  of  the  palm  are 
also  given  in  Table  1. 
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Figure  10.  Complex  permittivity  of  the 
skin  of  the  palm. 
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Effect  of  Millimeter-Wave  Irradiation  on  Growth  of 

Saccharomyces  cerevisiae 

LUCIANO  FURIA.  member,  ieee.  DOUGLAS  W.  HILL,  and  OM  P  GANDHI,  fellow,  ieee 


Abstract— Cultures  at  Saccharomyces  ctrtnsiat  were  exposed  for  4  h 
to  millimeter  waves  in  three  frequency  ranges  between  41.650  and 
41.798  GHz.  The  irradiation  frequency  was  stabilized  to  within  — 50 
Hz.  The  temperature  difference  between  irradiated  and  sham-irradi¬ 
ated  samples  was  maintained  to  within  ;0.01"C.  Growth  was  mea¬ 
sured  optically  during  the  irradiation,  and  viability  counts  were  done 
at  the  end  of  the  irradiation.  At  least  three  experiments  were  per¬ 
formed  at  each  of  15  frequencies. 

No  differences  larger  than  ;4  percent  were  detected  in  the  growth 
rates  at  any  of  the  selected  frequencies.  Such  differences  were  not  sig¬ 
nificant  at  the  95  percent  confidence  limit.  Results  obtained  with  plate 
counts  correlated  favorably  with  the  optical  absorbance  data.  While 
our  data  are  in  contrast  with  those  reported  from  other  investigators, 
these  experiments  support  conclusions  of  our  previous  studies,  and  of 
some  other  investigators,  showing  that,  under  strictly  controlled  con¬ 
ditions.  no  statistically  significant  nonthermal  effects  can  be  induced 
by  millimeter-wave  irradiation  uf  a  variety  of  prokaryotic  and  eukar¬ 
yotic  cells. 

Introduction 

SINCE  1968  several  investigators  have  reported  a  va¬ 
riety  of  biological  effects  induced  by  irradiation  with 
millimeter  waves  (mm-waves).  Many  have  further 
claimed  strong  dependence  on  the  irradiation  frequency. 
Representatives  of  these  findings  include  effects  on  growth 
rates  of  Rhodotorula  rubra  [  1  ],  Saccharomyces  cerevisiae 
|2)-|6|  Escherichia  coli  [7].  [8|.  Candida  albicans  { 1  ].  (9). 
reduced  viability  in  Sacch.  cerevisiae  (10|.  protection  of 
rabbit  bone  marrow  cells  from  X-ray  damage  (11|.  in¬ 
creased  coltcin  ( 12]  and  lambda  phage  induction  in  £.  coli 
[131.  and  puffing  of  giant  chromosomes  of  Acricotopus 
htadus.  an  insect  [  14J. 

Some  repons  indicated  changes  in  the  magnitude  of  the 
etfect  when  the  frequency  was  shifted  by  a  few  mega¬ 
hertz.  Furthermore,  some  reports  were  characterized  by 
the  low  power  densities  needed  to  induce  the  effect:  a 
power  density  as  low  as  10  pW/cm;  was  reported  to  effect 
colicin  induction  in  £.  coli  (12).  An  extensive  review  of 
the  published  literature  can  be  found  m  [15]. 

Unfortunately,  many  of  the  above  mentioned  reports 
lack  essential  details  of  the  experimental  procedure  and 
data,  making  independent  duplication  or  evaluation  quite 
difficult  Some  of  tnem.  however,  have  been  indepen- 
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dently  repeated,  mostly  with  negative  findings.  Some  o 
the  negative  experimental  results  include:  no  effects  ot 
growth  rate  of  £.  coli  [16]  and  no  changes  in  viability  o 
Sacch.  cerevisiae  [17],  No  evidence  of  cvtological  effect: 
as  determined  by  electron  microscopy  [18).  and  in  proteii 
synthesis  [19]  were  observed  in  BHK-21  cells  exposed  t< 
202  frequencies  in  the  £-  and  II-  band  at  various  powe: 
densities.  Finally,  no  statistically  significant  effects  wen 
found  on  colicin  induction  [20],  [21J.  on  mutation  rate: 
of  Salmonella  ryphimurium  or  on  induction  of  lambdi 
phage  in  £.  coli  bv  U-  and  £-band  mm- wave  irradiatioi 
[22|.  [23], 

A  strong  dependence  on  irradiation  frequency  of  thi 
growth  rate  of  yeast  suspension  has  been  reported  mde 
pendentiy  by  two  research  groups.  Devvatkov  [1|  re 
ported  increases  up  to  +30  percent  in  the  growth  rate  o 
R.  rubra  exposed  for  15  h  at  7. 18  mm  (41.783  GHz)  an< 
a  decreased  growth  rate  of  up  to  -40  percent  when  ex 
posed  at  7.19  mm  (41.725  GHz).  However,  no  indica 
tions  were  given  about  expenmental  procedures  such  a: 
irradiation  conditions,  handling  of  the  controls,  tempera 
ture  control  and  measurement,  frequency  stability,  am 
power  density. 

Grundler.  Keilmann.  and  co-workers  [2],  (b|  have  aKo 
reported  that  the  irradiation  of  Sacch.  cerevisiae  in  thi 
4 1.650-4 1.800  GHz  range,  showed  both  growth  enhance¬ 
ment  (up  to  +15  percent)  and  growth  inhibition  idown  to 
-13  percent)  frequency  regions.  Frequencies  exhibiting 
significantly  different  effects  were  reported  to  be  only  j 
MHz  apart.  The  data  suggest  that  the  yeast  suspensions 
showed  a  sensitivity  to  frequency  changes  of  less  than  I 
part  in  5000.  Due  to  this  reported  frequency  dependence, 
this  phenomenon  has  been  termed  a  "resonant  effect”  1 21- 
[6]. 

Although  such  reported  frequency  sensitivity  has  stim¬ 
ulated  interest  in  the  possible  mechanism  and  suet  si  o' 
action  of  mm-waves  |I8|-|23|.  !24|.  (25).  it  is  essentia 
to  establish  the  validity  of  these  observations  and  then 
interpretation  In  the  following  we  report  on  our  experi¬ 
ments  designed  to  study  further  the  possible  effects  ot  mm- 
waves  on  the  growth  ot  Sacih  icrcvmac 
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glucose  agar  i Neopeptone*  10  g/l.  glucose  40  g  1.  Bacto 
agar*  15  g  1)  stored  at4  =  C  and  passed  even  month.  From 
these  cultures,  organisms  were  inoculated  in  Sabouraud 
glucose  broth  (Neopeptone*  10  g«l.  glucose  20  g/l.  pH 
5.6)  and  incubated  overnight  at  32°C  on  an  orbital  shaker 
(40-50  rpm).  The  overnight  culture  was  diluted  1:400 
and  incubated  at  32 ®C  for  2  h  on  an  orbital  shaker  before 
loading  each  of  two  iiradiation  circuits  with  2.5  ml  of  the 
suspension.  The  irradiation  and  sham  irradiation  were 
started  after  another  hour  of  incubation  in  the  circuits  so 
as  to  irradiate  the  yeast  during  the  log  phase  of  growth. 
The  temperature  difference  between  irradiated  and  sham- 
irradiated  suspensions  was  maintained  within  ±0.01  °C 
using  a  nichrome  heater  placed  under  the  sham-irradiated 
sample  holder.  The  suspensions  were  recirculated  with  a 
peristaltic  pump  LKB  2115  with  a  flow  rate  of  7  ml/min. 
The  experimental  temperature  was  32.0  ±  0.1  °C.  The 
temperature  coefficient  of  the  growth  rate  was  measured 
between  30  and  34 °C  and  was  0.024  h_l/°C.  The  mean 
growth  rate  at  32°C  was  0.595  h"1.  corresponding  to  a 
duplication  time  of  70  min. 

The  irradiation  time  was  4  h  with  power  in  the  wave¬ 
guide  set  at  20.0  ±  0.5  mW.  Due  to  the  tapered  design 
of  the  sample  holder,  the  power  coupling  of  mm-wave 
radiation  to  the  suspension  was  larger  than  99  percent  as 
estimated  from  reflected  power  measurements.  Even 
though  we  used  continuous  wave  (CW)  microwaves,  cir¬ 
culating  the  suspension  did.  in  facl.  provide  a  modulation 
of  sorts  with  a  period  of  30  s  and  a  duty  cycle  of  0.06. 

Small  flow-through  cuvettes  were  placed  within  the  re¬ 
circulating  circuits  to  allow  turbidity  measurements.  For 
both  the  irradiated  and  sham-irradiated  samples,  triplicate 
->et.s  of  data  were  collected  at  the  beginning  of  the  irradia¬ 
tion  and  hourly  thereafter  All  the  data  collection,  stor¬ 
age.  and  analysis  was  performed  with  an  IBM  personal 
computer.  Absorbance  values  at  the  observation  time  in¬ 
tervals  were  analyzed  with  a  regression  analysis  program 
to  tit  an  exponential  curve.  The  output  of  the  program 
gave  the  estimated  mean  values  of  growth  rates  n,  and  n,. 
tor  irradiated  and  sham- irradiated  suspensions.  The  pro¬ 
gram  also  indicated  the  95  percent  confidence  limits  of 
the  estimated  growth  rates  for  each  experiment.  A  typical 
exponential  growth  curve,  measured  optically,  for  an  ex¬ 
periment  without  either  mm-wave  radiation  or  compen¬ 
satory  electric  heating  is  shown  in  Fig.  1. 

Plate  counts  were  done  at  the  end  of  the  irradiation  ex¬ 
periment  using  three  0  2  ml  aliquots  from  each  suspen¬ 
sion  After  tenfold  serial  dilution,  the  suspensions  were 
plated  in  quadruplicate  and  colonies  counted  after  48  h  of 
incubation  at  32°C  Dilutions  were  chosen  to  give  ap¬ 
proximately  200  colonies  per  plate 

B  Microwave  Svstem 

An  mentioned  previously,  the  main  feature  of  a  so- 
called  "resonance  effect"  is  its  Ntrong  dependence  on  the 
irradiation  frequency  Therefore.  Mime  of  the  basic  re¬ 
quirement  tor  ob\erv  mg  nunA  an  effect  Nhould  be  I  i  high 
frequency  Ntabtlitv  within  a  single  experiment.  2)  high 
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Fig.  1  A  typical  growth  curve  of  Sacchuromvcrs  cerevisi nr  suspension 
measured  optically  in  the  absence  of  mm-waves  For  this  particular  ex- 
penment  the  growth  rate  of  the  "sham-irradiation  '  was  0.53 1  h  .  while 
for  the  "irradiauon"'  circuit  tt  was  0522  h  1  The  95  percent  ..onridence 
limits  were  about  ±0.015  h"’  (±2.5  percent)  The  points  represent  the 
average  of  three  measurements  taken  at  the  indicated  time  intervals,  me 
lines  represent  the  titled  curves  calculated  with  an  exponential  regression 
program.  The  cell  concentrations  at  the  start  and  at  the  end  of  an  exper¬ 
iment  were  in  the  6.5  x  I05- I  I  x  10“  and  T  x  10"- 1  2  x  io’  cell  ml 
ranges,  respectively  The  time  intervals  and  dilution  lactorv  were  chosen 
so  as  to  do  experiments  only  dunng  the  log  phase  ot  ihe  > east  growth 
and  to  use  the  linear  response  region  of  the  photometric  svsiem 

precision  in  the  frequency  resetting  from  experiment  to 
experiment,  3)  high  accuracy  in  the  absolute  frequency 
determination,  and  4)  low  noise  and  narrow  bandwidth  of 
the  source.  After  testing  various  mm-wave  sources  such 
as  free  running  klystrons  and  IMPATT  diodes,  we  chose 
a  klystron.  Varian  model  VA302-BT.  phase-locked  by  a 
source-locking  microwave  counter.  EIP  model  578  By 
using  a  highly  stable  quartz  crystal  enclosed  in  a  temper¬ 
ature  controlled  oven,  the  rated  frequency  stability  of  the 
counter  was  ±50  Hz  for  a  center  frequency  around  40 
GHz.  This  stability  was  for  at  least  4  h.  the  typical  length 
of  an  experiment.  The  absolute  frequency  could  be  mea¬ 
sured  with  up  to  11  significant  digits.  The  rated  aging  of 
the  reference  quartz  crystal  is  claimed  bv  the  manufac¬ 
turer  to  give  a  long-term  accuracy  ( 1  year)  of  the  absolute 
frequency  of  less  than  ±8  kHz.  The  residual  half-power 
bandwidth  was  40  kHz  when  phase-locked,  as  measured 
with  a  Tektronix  491  spectrum  analyzer. 

The  CW  power  in  the  waveguide  was  measured  with  a 
thermistor.  Hughes  model  4489H.  connected  to  a  power 
meter,  Hewlett-Packard  model  432A.  A  variable  atten¬ 
uator  was  inserted  before  the  thermistor  to  cut  off.  from 
time  to  time,  the  microwave  power  to  the  thermistor, 
therefore  allowing  correction  for  possible  thermal  drift  of 
the  thermistor.  The  sample  holders  were  mounted  at  a  6 3 
angle  directly  on  the  waveguides  i  Fig.  2)  thereby  provid¬ 
ing  enough  tapering  to  ensure  a  measured  power  coupling 
of  jbout  99  percent.  It  was  confirmed  during  preliminary 
tests  that  the  reflected  power  was  always  less  than  1  per¬ 
cent.  regardless  of  the  irradiation  frequency  The  fixed 
geometry  of  the  sample  holder  did  not  require  routine 
measurement  ot  the  reflected  power.  A  complete  block 
diagram  of  the  microwave  circuit  is  shown  in  Fig.  3 
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Fig.  2.  Sample  holder  used  tor  irradiation.  The  Terton-*  coated  thermo¬ 
couples  were  inserted  in  the  suspensions  through  holes  in  the  silastic 
rubber  tubing.  A— WR-19  waveguide:  fl— cell  suspension  circuit:  C— 
water  jacket  circuit:  0—  thermocouple  for  the  temperature  measurement 
of  the  sample  at  the  exit  of  the  sample  holder.  Another  thermocouple 
was  placed  also  at  the  sample  holder  inlet.  The  directions  of  how  are 
indicated  by  the  arrows  on  the  circuits.  The  nichrome  heater  was  placed 
between  the  waveguide  and  the  sham-irradiation  sample  holder. 


Fig  3  Block  diagram  ot  the  microwave  system.  The  section  necessary  for 
the  reflected  power  measurement  was  removed  during  the  actual  exper¬ 
iments. 


C.  The  Irradiation  Chamber 

The  chamber  was  made  with  an  aluminum  box  lined 
with  polyurethane  foam.  On  the  bottom  of  the  chamber 
we  placed  a  copper  tank  in  which  temperature  controlled 
water  was  recirculated  at  the  rate  of  10  1/min.  The  water 
was  circulated  from  a  water  bath  kept  at  32.0  ±  0.1  °C 
by  a  YSI  72  proportional  temperature  controller  and  a 
solid-state  cooler  (Whirlpool).  The  waveguides  with  the 
sample  holders  were  bolted  on  the  copper  tank.  To  obtain 
a  uniform  temperature  distribution  we  found  it  essential 
to  place  an  electric  fan  inside  the  chamber.  The  60  Hz 
magnetic  field  generated  by  this  fan  did  not  differ  from 
the  background  intensity.  0  03  G.  at  either  the  sample 
holder  or  the  cuvette  location  A  top  view  of  the  irradia¬ 
tion  chamber  is  shown  in  Fig  4  The  ulass  sample  holder 
nad  a  volume  of  i)  25  ml  '2  5  -  0  5  <02  cmi  and  was 
connected  to  the  peristaltic  pump  with  silastic  rubber  tub¬ 


Fig.  4.  Top  view  of-the  irradiation  chamber.  The  copper  tank  is  at  thi 
bottom  of  the  chamber.  A  metallic  lid  was  placed  on  the  chamber  dunni 
the  experiments  A— peristaltic  pump:  B— electnc  fan.  C— cuvettes  fa 
turbidity  measurements:  D— optical  liber  bundles:  £— sham-irradiatiot 
sample  holder.  F—  irradiation  sample  holder:  and  G—  temperature  prow 
for  water  jackets.  I 


ing  (1.3  mm  i.d.).  Water  jackets  were  placed  around  thi 
glass  sample  holders  to  improve  the  thermal  exchange 
Water  for  the  water  jackets  was  circulated  at  a  rate  of  1  I 
min,  from  the  same  water  bath  supplying  the  copper  tank 
Teflon*  coated  thermocouples  (Bailey-Sensortek  IT- 18 
were  placed  in  the  suspensions  at  the  inlets  and  outlets  c 
the  sample  holders.  These  were  connected  to  a  digital  dil 
ferential  thermometer  (Bailey-Sensortek  TH-6D)  that  a 
lowed  a  rated  precision  in  the  absolute  mode  of  ±0. 1  °( 
and  of  ±0.01  °C  in  the  differential  mode.  The  temperatur 
difference  was  minimized  by  adjusting  the  current  in  th 
nichrome  heater.  No  temperature  increase  due  to  the  m; 
crowave  or  electrical  heating  was  detected  in  the  susper 
sions  throughout  the  experiment.  The  use  of  small  ci 
vettes.  HELLMA  144  O- NS.  and  the  placement  of  the  hea 
of  the  peristaltic  pump  inside  the  irradiation  chamber  a 
lowed  the  use  of  a  total  suspension  volume  of  2.5  ml. 

D.  Optical  System 

A  fiber  optic  system  was  used  to  measure  the  turbidit 
of  the  yeast  cultures  and  in  turn  the  growth  rates  of  th 
organisms.  We  used  optical  fibers  made  with  Croton 
(Dupont)  fiber  bundles  (48  fibers,  each  0.25  mm  in  d 
ameter).  An  interference  filter  (Edmund  Scientific  Coi 
poration.  650  nm  center  wavelength.  10  nm  bandwidth 
was  placed  between  the  optical  fibers  and  a  150  W  froste 
tungsten  lamp.  A  stabilized  ac  power  supply  was  used  I 
obtain.a  light  output  stable  to  within  ±0.5  percent.  Th 
turbidity  was  determined  with  the  optical  fibers  placed  i 
contact  with  the  walls  of  the  cuvettes.  The  incident  ligl 
intensity  at  the  cuvette  was  about  6.5  jxW.  while  the  tram 
mitted  light  intensity  ranged  between  330  and  230  nM 
The  light  intensity  was  measured  with  a  Photodvne  XU 
photometer/ radiometer  with  sensor  head  model  250. 

Rest  lts 

The  possible  effects  of  mm-waves  on  yeast  suspensioi 
were  monitored  by  measuring  the  growth  rate  ot  the  cvi 
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RELATIVE  CONCENTRATION 

Fig  5  A  typical  plot  ot  absorbance  versus  concentration.  A  concentranon 
ol  I  corresponds  to  a  1  10  dilution  of  the  overnight  culture  The  points 
represent  the  average  ol  three  measurements  taken  at  the  indicated  con¬ 
centrations.  The  slope  for  the  "sham-irradiation-'  circuit  had  a  coeffi¬ 
cient  of  0.112  with  a  45  percent  confidence  limit  of  ±0  003  The  "ir¬ 
radiation''  circuit  slope  had  a  coefficient  of  0.113  f  0.002.  The  lines 
represent  the  fitted  curves  calculated  with  a  linear  regression  program. 
The  shaded  bar  indicates  the  absorbance/  concent  rat  ion  range  usually  used 
for  the  experiments. 


tures  and  by  performing  a  plate  count  at  the  end  of  each 
of  the  experiments.  The  absorbance  was  calculated  as  Abs 
=  log  ( PJP„ )  where  Pm  was  the  light  intensity  transmit¬ 
ted  through  the  reference  cuvette  filled  with  Sabouraud 
glucose  broth  and  P,r  was  the  light  intensity  transmitted 
through  the  cuvettes  filled  with  either  the  sham-irradiated 
or  the  irradiated  suspension.  Before  each  experiment,  the 
light  transmitted  by  the  three  cuvettes,  all  filled  with  the 
growth  medium,  was  verified  to  be  within  the  resolution 
of  the  photometer  ii  e..  I  count  =  1  nW|. 

The  linearity  of  the  response  of  the  optical  system  of 
both  circuits  was  determined  before  each  experiment  using 
a  twofold  dilution  senes,  whose  cell  concentration  range 
was  close  to  that  of  the  actual  expenment.  A  typical  ab¬ 
sorbance  versus  concentration  curve  is  shown  in  Fig.  5. 
For  a  total  of  38  tests,  the  average  percentage  difference 
between  the  growth  rates  measured  in  the  two  circuits  was 
-0.77  percent,  with  a  standard  deviation  of  0.57  percent. 
The  maximum  difference  ever  detected  was  +2.5  per¬ 
cent. 

The  growth  rate  was  obtained  by  using  a  regression 
analysis  to  fit  the  data  to  an  exponential  curve  C(r)  - 
C/'.  where  C(r)  is  the  concentration  at  time  r.  C„  is  the 
initial  concentration,  and  m  is  the  growth  rate  in  h-'  The 
percentage  difference  in  the  growth  rate  of  veast  cells  was 
calculated  as  100  *  i  g  ~  n,  )'m.  ■  where  n  and  n,  are 
the  growth  rates  for  irradiated  and  sham-irradiated  cul¬ 
tures.  respectively 

We  performed  lb  experiments  without  either  mm-wave 
irradiation  or  electrical  heating,  and  under  such  condi¬ 
tions  the  average  difference  in  growth  rates  between  the 

irradiated "  and  the  "■'ham-irradiated"  circuit  was  -  I  5 
r  3  3  percent  The  average  ^5  percent  confidence  limit 
u  the  estimated  arowth  rule  was  about  -2  3  percent  re- 
cjrdless  of  the  circuit  or  ot  the  presence  ot  mm-wjves 
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Fig.  6.  Effects  of  millimeter-wave  irradiation  on  the  absolute  growth  rate 
of  Saccharnmyces  cerevisiae  as  a  function  of  frequency  lor  an  absorbed 
power  of  20  ±  0.J  mW  Each  point  and  bar  represent  the  average  -  1 
SO  of  the  experiments  at  each  frequency  Unless  otherwise  indicated  tn 
parentheses,  three  experiments  were  performed  at  each  frequency 


From  these  data  we  can  infer  that,  with  at  least  three  ex¬ 
periments  for  each  frequency,  the  minimum  difference  L. 
significant  at  the  95  percent  confidence  limit,  was  about 
±3.5  percent  as  calculated  from  L  =  1.96  u/vi i.  where 
a  is  the  average  standard  deviation.  3.2  percent,  and  n  is 
the  number  of  experiments  [26].  In  the  three  frequency 
ranges  that  were  covered,  the  selected  frequencies  were  2 
MHz  apart.  At  least  three  experiments  were  performed  at 
each  frequency;  the  actual  numbers  are  indicated  in  Figs. 
6  and  7. 

In  Fig.  6  the  results  of  mm-wave  irradiation  on  the 
growth  rate  of  yeast  measured  optically  are  shown  as 
averages  of  the  absolute  growth  rates  n  and  u.  •  In  Fig. 
7(a)  the  same  data  are  also  plotted  a*  a  percentage  differ¬ 
ence  of  the  growth  rates.  Paired  i  tests  were  performed  on 
the  absolute  growth  rates  of  irradiated  versus  sham-trra- 
dialed  cultures,  grouped  by  frequency  No  statisticallv 
significant  differences  were  detected  ( P  >  0.05)  at  any  of 
the  selected  frequencies. 

For  most  of  the  experiments  we  also  performed  a  plate 
count  at  the  end  of  irradiation.  The  percentage  differences 
between  the  number  of  viable  population  units  of  the  ir¬ 
radiated  culture  minus  those  of  the  sham-irradiated  are 
shown  in  Fig.  7(b).  No  statistically  significant  differences 
(P  >  0.05)  were  detected  by  unpaired  t  test  statistics  per¬ 
formed  for  each  experiment,  except  for  some  experiments 
at  some  of  the  frequencies  at  which  0.01  <  P  <  0  05 
This,  however,  never  occurred  for  more  than  half  of  the 
total  number  of  experiments  done  at  each  frequency 

We  have  also  performed  some  experiments  at  higher 
power  levels  at  a  fixed  frequency  of  41  698  GHz  The 
results  are  showed  in  Fig  8  No  effects  were  detected 
even  at  these  higher  powers. 

Discussion 

In  the  past  some  investigators  indicated  that  the  mm- 
wave  1 1 1  at  different  trequenev  ranges  alters  the  crowth 
rates  ot  R.  rubra  |l|.  Saa  li.  i  crcvtsiac  |2|-|b|.  E.  >  <>u 
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Fig.  7.  Effects  of  millimeter-wave  irradiation  on  Saccharomyces  cerevis- 
iae  as  a  function  of  frequency,  (a)  Percentage  differences  ol  the  growth 
rates,  (b)  Percentage  difference  in  the  number  of  viable  populations  units. 
Each  point  represents  the  average  ±1  SD  of  the  experiments  at  each 
frequency  Unless  otherwise  indicated  in  parentheses,  three  experiments 
were  performed  at  each  frequency. 
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Fig  8  Effects  of  higher  power  levels  on  ihc  growth  rate  ot  yeast.  Fre- 
4uenc\  rixed  at  41  6980  GHz.  The  points  at  40  and  80  mW  represent 
the  result  of  single  experiments 

[7],  and  C.  albicans  [9],  Other  reports  also  indicated  that 
mm- waves  enhance  the  induction  of  colicin  [12]  and  of 
lambda  phage  [13]  in  E.  coli.  These  reports  prompted 
several  comparative  studies  on  prokaryotic  and  eukaryotic 
cells  such  as  E.  coli  [16],  Sacch.  cerevisiae  [17].  BHK- 
21  [18].  [19],  and  5.  typhimurium  [22],  [23].  All  these 
studies  failed  to  indicate  any  nonthermal  effect  of  mm- 
waves.  In  particular,  a  study  on  protein  synthesis  in  mam¬ 
malian  cells  BHK-21  exposed  at  various  power  levels 
failed  to  reveal  any  change  due  to  mm-waves  at  202  fre¬ 
quencies  in  the  38-48  GHz  and  65-75  GHz  ranges  [19], 
As  shown  in  Figs.  6  and  7  the  present  study  supports 
the  conclusions  that  mm-waves  do  not  affect  in  a  detect¬ 
able  manner  the  growth  rate  or  the  viability  of  yeast  cells. 
The  average  growth  rate  difference  of  irradiated  cultures 
of  Sacch.  cerevisiae  differed  by  no  more  than  ±4  percent 
from  the  sham-irradiated  at  any  of  the  selected  frequen¬ 
cies,  while  the  plate  counts  did  not  differ  by  more  than 
*  15  percent.  The  larger  differences  in  the  plate  count  may 
nave  been  due  to  flight  initial  concentration  differences 
resulting  from  sampling  errors  that  might  have  occurred 
when  the  circuits  were  loaded  at  the  beginning  of  the  ex¬ 


periment.  This  difference,  however,  would  have  not  re¬ 
sulted  in  significant  differences  in  the  growth  rates. 

The  spectral  characteristics  of  our  microwave  source 
can  be  favorably  compared  to  those  of  a  similar  experi¬ 
ment  repotted  by  Grundler  et  al.  [6],  where  the  bandwidth 
was  ±0.5  MHz.  the  frequency  stability  was  ±0.3  MHz. 
and  the  absolute  frequency  was  accurate  to  within  ~0. 1 
MHz.  In  particular,  the  use  of  very  accurate  frequency 
counters  in  both  this  study  and  the  ones  reported  in  [5], 
[6]  excludes  the  possibility  that  we  might  have  used  ab¬ 
solute  frequencies  offest  by  more  than  a  few  hundreds  of 
a  kilohertz.  Since  a  systematic  scanning  of  even  a  small 
portion  of  the  millimeter-wave  band.  30-100  GHz.  with 
a  frequency  step  of  2  MHz  would  have  required  an  inor¬ 
dinate  amount  of  time,  we  have  focused  the  experiments 
at  the  frequency  ranges  in  which  the  largest  effects  have 
been  reported  in  the  past  on  the  growth  of  Sacch.  cere¬ 
visiae  [2]— [6] . 

The  objection  could  be  made  that  the  difference  in  the 
results  could  be  ascribed  to  the  differences  in  strain.  This 
objection  postulates  the  existence  of  a  wild-type  strain  of 
yeast  that  is  naturally  sensitive  to  mm-wave  radiation.  On 
the  basis  of  experience  with  induced  sensitivity  to  physi¬ 
cal  agents,  like  UV  for  example,  it  seems  extremely  un¬ 
likely  that  this  could  happen.  In  the  case  of  UV,  for  ex¬ 
ample,  the  resistance  is  obtained  by  selecting  those  cells 
(usually  1  in  106)  that  survive  UV  treatment.  In  our  case 
the  existence  of  a  mm-wave  sensitive  strain  would  imply 
that  either  a)  the  strain  used  by  Grundler  and  Keilmann 
has  been  used  after  becoming  sensitive  to  exposure  to  mm- 
wave  radiation,  or  hi  we  used  a  strain  that  tor  some  reason 
had  lost  its  sensitivity  to  mm-wave  radiation.  The  prob¬ 
ability  of  either  phenomenon  happening,  even  though  the 
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exact  probability  cannot  be  determined,  seems  to  be  ex¬ 
ceedingly  small.  Furthermore,  biological  effects  of  mm- 
waves  have  been  reported  in  a  variety  of  different  living 
organisms. 

It  should  be  pointed  out  that  the  points  and  the  bars 
shown  by  Grundier  et  ai.  in  Fig.  9  of  [5]  do  not  represent 
the  average  and  the  standard  error,  respectively,  of  mul¬ 
tiple  experiments  at  the  same  frequency.  Instead,  the 
points  are  the  results  of  single  experiments  and  the  “error 
bar”  was  obtained  as  follows.  Two  different  spectropho¬ 
tometers  monitored  the  growth  of  two  cultures  prepared 
from  the  same  suspension.  The  maximum  percent  differ¬ 
ence.  over  several  experiments  without  microwaves,  be¬ 
tween  the  growth  rates  of  the  two  cultures  was  then  used 
'as  the  estimated  error.  Therefore,  the  "error  bars”  in  Fig. 
9  of  [5]  all  duplicate  the  same  information.  This  repre¬ 
sentation  also  implies  that  the  authors  of  [5],  [6]  assume 
the  estimated  experimental  error  to  remain  constant  re¬ 
gardless  of  the  irradiation  frequency.  If  all  the  points 
within  ±4  percent,  i.e.,  their  estimated  error,  of  Fig.  9 
in  [51  are  considered  as  the  experimental  “noise."  then 
the  use  of  a  three-point  interpolated  curve  to  connect  all 
the  experimental  data  with  a  resonant-like  curve  becomes 
quite  unjustified.  However,  even  though  in  (5).  [6]  there 
are  not  enough  experimental  data  at  most  of  the  selected 
frequencies  to  claim  a  definite  effect  on  the  growth  rate, 
there  is  an  apparent  similarity  between  the  results  ob¬ 
tained  with  two  independent  sets  of  experiments,  using 
two  different  irradiation  systems.  Furthermore,  this  sim¬ 
ilarity  is  also  supported  by  the  cross  correlation  analysis 
of  the  two  sets  of  experiments  [5],  [6], 

The  experiments  reported  in  [2]-(6]  were  not  done  at 
constant  mm-wave  power,  instead,  the  power  varied  from 
experiment  to  experiment  between  5  and  40  mW.  intro¬ 
ducing  thereby  another  uncontrolled  variable  in  the  ex¬ 
perimental  procedure.  However,  both  in  the  present  study 
and  in  the  ones  reported  earlier  [2]-[6],  no  cell  was  ex¬ 
posed  continuously  to  mm-wave  energy  due  to  the  small 
skin  depth  of  mm-waves  and  either  recirculation  or  stir¬ 
ring  of  the  cell  suspension.  In  fact,  two  of  the  more  not¬ 
able  differences  between  the  study  in  (2]-[6j  and  the  pres¬ 
ent  one  are  a)  the  average  time  for  a  single  exposure  (0.06 
s  versus  2.0  s  for  the  present  study);  b>  the  average  inter¬ 
val  between  subsequent  exposures  (1.2  s  versus  30  s). 
Nevertheless,  the  total  time  that  the  average  cell  was  ex¬ 
posed  to  mm-wave  energy  is  comparable  (5  percent  ver¬ 
sus  6.7  percent  for  the  current  study). 

As  a  suggestion  for  future  studies  of  biological  effects 
of  mm-waves.  we  would  like  to  point  out  that  in  the  past 
when  the  power  used  was  less  than  5  mW/cm:  such  effects 
have  been  considered  "nonthermar '  in  nature  [2|-|8), 
(I3|  [I7].  suggesting  that  the  microwave-induced  heat¬ 
ing  could  not  have  caused  the  reported  effect.  However, 
the  skin  depth  6  is  very  small  in  lossy  materials  such  as 
biological  media  containing  large  amounts  of  water  For 
example  in  distilled  water,  o  is  between  0.7S  and  0  23 
mm  tor  rreuuencies  m  the  range  50-500  GHz  This  leads 
to  xery  nign  '•.dues  of  the  specific  absorption  rate  (SAR) 
in  the  surface  layers  of  the  exposed  medium  even  tor  ap¬ 


parently  small  incident  power  densities.  For  example,  at 
I  raW  cm',  the  surface  SAR  in  water  at  40  GHz  is  13  - 
W  kg.  On  the  other  hand,  the  power  density  P  decreases 
exponentially  as  P  ~  P  e~~  ~'  where  P  is  the  inciuent 
power  density  and  r  is  the  distance  from  the  surface  ex¬ 
posed  to  the  incident  field.  Therefore,  any  organism  lo¬ 
cated  more  than  3-4  skin  depths  from  the  exposed  surface 
is  subjected  only  to  a  negligible  fraction  of  the  incident 
mm-wave  radiation. 

It  is  therefore  important  to  determine,  on  an  individual 
basis,  whether  such  high  values  of  power  deposition  can 
lead  to  subtle  thermal  effects  or  not.  Some  of  the  reported 
effects  have  been  more  appropriately  called  nonthermai 
when  either  of  the  following  circumstancs  occurred;  ai 
the  microwave  irradiation  induced  an  effect  opposite  to 
that  expected  by  a  comparable  increase  in  temperature:  bi 
when  the  effect  appeared  to  be  strongly  dependent  on  the 
irradiation  frequency. 

Conclusions 

The  current  experiments  were  specifically  designed  to 
gather  data  that  would  help  establish  the  presence  or  ab¬ 
sence  of  nonthermai  effects  and  frequency  specific  effects 
of  mm-waves  on  the  growth  of  cells.  This  effect,  first  re¬ 
ported  in  1968  [7]  for  the  growth  of  E.  coli  was  followed 
by  other  reports  indicating  specific  induction  of  lysogenic 
E.  coli  [13],  colicin  induction  [12],  and  others  [lJ-[6], 
[8],  [9].  The  fundamental  and  empirical  significance  of 
these  alleged  frequency  specific  actions  of  nonionizing  ra¬ 
diation  in  the  mm-wave  range  hardly  needs  any  emphasis. 

Contrary  to  the  previous  reports  [21— [6]  from  other  in¬ 
vestigators.  millimeter  waves  did  not  induce  any  detect¬ 
able  effect  on  either  the  growth  rate  or  the  viability  of 
yeast  cells  exposed  for  4  h  to  ultrastable  millimeter  waxes 
between  41.650  and  41  "98  GHz.  This  report  adds  to  a 
growing  list  of  others  [151-[23]  showing  that  under  care¬ 
fully  controlled  experimental  conditions,  no  nontheimai 
effects  of  millimeter  waves  on  unicellular  organisms  are 
evident. 

This  conclusion  suggests  extreme  caution  when  uncon¬ 
firmed  reports  are  used  as  a  basis  for  broad  generaliza¬ 
tions  of  biological  effects  of  millimeter  waves  [2-4).  [25] 
by  authors  who  also  suggest  the  possible  use  of  mm-waves 
for  cancer  diagnosis  [27],  [28]  and  therapy  [29),  and  a 
possible  means  of  cell  communication  [24], 
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Cultures  of  Bacillus  megaterium  have  been  examined  with  laser -Ram an  spectroscopy  in  the  20-300  and  800-1800  cm'1 
ranges.  No  lines  were  detected  in  the  20-300  cm'1  range,  and  evidence  is  presented  that  sporadic  lines  observed  in  the 
other  range  can  be  ascribed  to  experimental  artifacts. 


The  presence  of  time -varying  lines  in  Raman  spectra 
of  cultures  of  Bacillus  megaterium  and  Escherichia 
coli  have  been  reported  [1—3]  in  recent  years.  Webb 
claimed  [3]  that  these  lines  could  be  observed  only  if 
the  cultures  were  in  an  active  metabolic  state  and  had 
been  “synchronized”  with  heat  and  cold  shocks  im¬ 
mediately  before  the  recording  of  Raman  spectra. 
Subsequently,  several  theoretical  reports  were  pub¬ 
lished  providing  possible  explanations  [4—6]  for  the 
reported  lines.  These  reports  suggested  either  a  relation 
between  phonon  condensation  and  cell  cycle  [4],  or 
the  presence  of  soliton  vibrations  in  biomolecules  [6] . 
or  just  an  empirical  relationship  among  the  reported 
lines  (5).  However,  a  recent  report  [7]  has  questioned 
their  presence  due  to  the  low  concentration  and  very 
low  Raman  scattering  cross  section  of  the  biomole¬ 
cules  possibly  involved  in  the  inelastic  scattering  pro¬ 
cess.  Furthermore,  in  [7],  experimental  evidence  was 
provided  that,  in  general,  an  inhomogeneous  bacterial 
culture  might  produce  large  fluctuations  in  Mie  scat¬ 
tering  during  the  first  minutes  after  resuspension,  and 
consequently  altering  the  baseline  of  the  recorded 
spectrum.  Webb’s  reports,  in  particular,  have  been 
used  by  Frohlich  [8,9]  and  Webb  himself  [3]  to  sup¬ 
port  Frohlich’s  hypothesis  [9-1 1]  concerning  the 
possible  presence  of  coherent  excitations  in  biological 
systems.  The  existence  of  these  coherent  oscillations 
has  been  conjectured  to  result  in  frequency-dependent 
biological  effects  induced  by  low-power  millimeter- 
wave  exposure  [9] ;  the  presence  of.  or  lack  of,  coher- 
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ent  oscillations  has  been  also  related  to  the  cancer  prob¬ 
lem  [8,12], 

We  report  Raman  experiments  with  cultures  of  B. 
megaterium  which  were  prepared  according  to  the  ex- 
perimental  procedure  reported  by  Webb  in  ref.  [3]. 
The  Raman  spectra  were  recorded  under  flowing  and 
nonflowing  conditions  either  at  room  temperature 
(21— 23°C)  or  at  a  constant  temperature  of  35°C.  It 
is  to  be  emphasized  that  in  refs.  [1-3]  no  mention  is 
made  about  any  kind  of  temperature  control  and  that 
the  spectra  were  recorded  under  nonflowing  conditions. 

Briefly,  the  B.  megaterium  cells  were  grown  for  48 
hours  in  nutrient  broth  at  3Z°C.  washed  three  times  in 
0.9  percent  NaCl  solution,  counted,  and  resuspended 
in  saline  at  a  concentration  of  1  X  107  cell/ml.  The 
cells  were  then  placed  in  a  water  bath  at  39°C  for  1 5 
min  followed  immediately  afterwards  by  an  ice  bath 
for  another  15  min.  The  suspension  was  then  centri¬ 
fuged  and  the  cells  resuspended  at  a  concentration  of 
2  X  107  ml-1  in  Davis  Minimal  Medium  (DMM) 
(K2HP04  7  g/2,  KH2P04  3  g/8.  Na3  citrate  0.5  g/8 
FeS04  0.01  g/8,  (NH4)2S04  1  g/8,  MgS04  +  7H20 
0.1  g/2,  glucose  2  g/8)  diluted  1 :3  in  double  distilled 
water  at  35°C.  All  the  suspending  media  were  previ¬ 
ously  filtered  through  a  0.22  m  Millipore®  filter.  For 
each  experimental  session,  a  single  bacterial  suspension 
of  about  10  ml  was  used  and  kept  in  a  water  bath  at 
either  35°C  or  at  room  temperature.  As  a  sample  hold¬ 
er.  we  used  a  NSG  Precision  cell  T73-FL  that  allowed 
us  to  circulate  both  the  sample  under  investigation 
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Fig.  1.  Low-wavenumber  Raman  spectra  of  B.  megaterium  cell 
suspension  after  a  heat  and  cold  shock  treatment  and  resus- 
penbon  in  nutrient  medium,  (a)  Spectrum  taken  under  flow¬ 
ing  condition  at  3S°C.  (b)  Spectrum  taken  at  room  tempera¬ 
ture  under  nonflowing  conditions.  The  optical  system  em¬ 
ployed  a  combination  of  a  double  monochromator  plus  a  third 
monochromator  with  a  spectral  bandwidth  of  4cm'1  and  with 
1  s  of  integration  time. 

and  temperature-controlled  water  in  an  outer  jacket. 
To  establish  sample  flow,  a  peristaltic  pump  set  at  S 
ml/ min  was  used. 

The  Raman  system  consisted  of  a  Spectra  Physics 
164/09  Argon-ion  laser  operating  at  4880  A.  a  Spex 
plasma  line  filter,  and  a  Spex  third  monochromator 
coupled  to  a  Spex  1403  double  monochromator.  The 
detection  system  was  composed  of  an  RCA  C3 1034 
cooled  GaAs  photomultiplier  and  a  Princeton  Applied 
Research  1 109  photon  counter.  The  whole  system  was 
controlled  by  an  LSI  1 1/23  minicomputer.  The  laser 
power  at  the  sample  was  about  150  mW  Typically, 
spectra  were  recorded  with  a  step  sue  of  4cm-1  and 
an  integration  tune  of  l  s. 

For  a  total  of  over  15  spectra,  no  lines  were  ever 
detected  in  the  20-300  cm-1  region  at  any  time  be¬ 
tween  10  mm  and  2  h  30  min  after  resuspension  in 
DMM.  at  either  room  temperature  or  at  35°C  and  un¬ 
der  flowing  or  nonflowing  conditions.  Typical  spectra 
are  presented  in  fig.  1 ,  For  over  50  spectra  in  the  800- 
1800  cm*1  range,  range,  we  frequently  observed 
random  lines  under  flowing  conditions  from  B  mega- 
terium  at  35°C  or  at  room  temperature  and  also  from 
the  suspension  medium  (DMM)  alone,  fig.  2.  Only 
rarely  did  we  observe  small  random  lines  from  non¬ 
flowing  preparations  Typical  spectra  under  nonlfow- 
ing  conditions  are  presented  in  fig.  3  Both  ranges  were 
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Fig.  1  Apparent  Raman  lines  shown  under  flowing  conditions 
in  spectra  of  either  (a)  B.  megaterium  ceils  suspended  in  nutri¬ 
ent  medium  (DMM  1:3)  at  35*C  after  heat  and  cold  shock,  (b) 
B.  megaterium  suspended  in  untnent  medium  after  heat  and 
cold  shock  at  room  temperature,  or  (c)  suspension  medium. 
Davis  minimal  medium  1:3.  The  optical  system  was  the  same 
as  for  spectra  of  fig.  1.  The  broad  peak  around  1640  cm-1  is 
due  to  the  H-0  bending  mode  of  water,  while  the  weak  peak 
at  980  cm*1  can  be  ascribed  to  the  presence  of  S04  group 
and  KjHPOa  and  KH3PO4  in  the  suspending  medium. 

also  examined  in  the  anti-Stokes  region  and  no  lines 
at  ail  were  ever  detected,  in  contrast  to  a  previous  re¬ 
port  [3). 

A  possible  explanation  for  these  lines  is  the  pres¬ 
ence  of  either  fluorescent  particles  or  cell  clumps  pass¬ 
ing  randomly  through  the  laser  beam.  As  these  random 
lines  can  be  reproduced  even  in  the  absence  of  biolog¬ 
ical  systems  of  any  kind,  we  thus  conclude  that  the 
observed  Raman  lines  do  not  rise  from  active  bacterial 
cultures.  It  is  suggested  that  experiments  be  performed 
on  simpler  and  more  charactenzable  biological  systems 
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Fig.  3  Raman  spectra  taken  under  nonllowmg  condition.  iai 
B  megaterium  cell  suspension  ibl  Davis  minimal  medium. 
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such  as  selected  antibody-antigene  complexes  and 
enzymatic  reactions.  With  these  systems,  m  fact,  it 
would  be  possible  to  use  greater  concentrations  and 
consequently  to  have  a  larger  total  scattering  signal. 
Furthermore,  the  lines  could  be  attributable,  at  least 
in  principle,  to  specific  biochemical  species,  thus  al¬ 
lowing  a  more  rigorous  test  of  any  hypothesis  regard¬ 
ing  possible  presence  of  coherent  excitations  in  biolo¬ 
gical  systems. 
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A  previous  paper  reporting  the  absence  of  biologically  related  lines  in  Raman  spectra  of  Bacillus  megatenum  has  been 
criticized.  A  reply  to  these  criticisms  is  presented. 


In  the  past  few  years  some  investigators  have  claim¬ 
ed  the  presence  of  lines  in  Raman  spectra  of  metabo¬ 
lizing  cells  while  no  lines  could  be  detected  in  spectra 
of  resting  cells  [1-3] .  Consequently  other  investiga¬ 
tors  have  linked  these  alleged  lines  to  internal  vibra¬ 
tions  of  Davydov  solitons  [4,5] ,  to  coherent  oscilla¬ 
tions  as  postulated  by  Frdhlich  [6]  and  to  nonlinear- 
optics  phenomena  [7,8] .  In  our  paper  [9]  we  present¬ 
ed  experimental  data  showing  that  no  biologically 
related  lines  were  found  in  Raman  spectra  of  Bacillus 
megatenum  suspensions  under  nonflowing  conditions 
in  the  20-300  and  800-1800  cm-1  regions  and  that 
under  flowing  conditions  artifacts  could  be  responsible 
for  random  lines  appearing  in  the  spectra  of  both  cell 
suspensions  and  medium  alone.  For  these  experiments 
we  closely  followed  the  experimental  protocol  given 
by  Webb  in  ref.  [3] . 

Del  Giudice  et  al.  [10]  have  raised  several  ques¬ 
tions  about  both  the  correctness  of  the  conclusions 
and  the  validity  of  experimental  data  presented  in  our 
previous  paper  [9] .  In  particular  our  conclusions  have 
been  attributed  more  to  an  alleged,  “...preconceived 
philosophy...”  and  to  “...motive  behind  the  two  re¬ 
cent  reports...”  rather  than  to  a  logical  connection  to 
the  experimental  data.  We  And  these  allegations  out 
of  place  in  a  scientific  forum,  especially  when  these 
are  totally  unsubstantiated.  After  a  very  careful  re¬ 
reading  of  our  paper  we  have  not  been  able  to  find 
any  suggestion  of  an  alleged  preconception.  Further¬ 
more  Del  Giudice  et  al.  attributed  to  us  a  criticism 
of  the  work  of  Webb  [1-3]  that  was,  instead,  point¬ 
ed  out  by  Cooper  and  Amer  [11].  Cooper  and  Amer 
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suggested  that  the  Raman  spectra  of  bacterial  ceils, 
prepared  according  to  Webb  [3]  were  likely  to  be 
featureless  due  to  the  relatively  low  cell  concentra¬ 
tions. 

In  the  present  reply  we  hope  to  offer  some  clarifi¬ 
cations  and  leave  the  conclusions  to  peers.  Some  of 
the  questions  asked  in  ref.  [10]  are,  however,  already 
explicitly  answered  in  our  original  paper  [9] .  In  the 
following,  citations  from  our  first  paper  [9]  are  in 
quotes: 

(a)  We  did  not  use  a  nonnutrient  medium.  The  leg¬ 
end  of  fig.  2  of  ref.  [9]  contains  a  printing  error:  “un- 
trient”  instead  of  nutrient.  We  really  wonder,  how¬ 
ever,  how  the  word  “untrient”  can  be  read  as  non¬ 
nutrient. 

(b)  The  experimental  conditions  of  35°C  and  room 
temperature,  were  used  because  the  latter  is  the  one 
used  by  Webb  in  ref.  [3]  and  the  former  is  an  op¬ 
timal  growth  temperature  for  B.  megaterium. 

(c)  In  addition  to  several  experiments  under  non¬ 
flowing  conditions,  we  used  also  a  flow-through  sys¬ 
tem,  since  this  was  also  the  recommendation  of  Drissler 
[12]  to  avoid  excessive  laser  heating  of  the  cells.  We 
pointed  out  in  our  paper  [9]  that  in  the  experiments 
reported  by  Webb  [1-3]  “...the  spectra  were  record¬ 
ed  under  nonflowing  conditions”. 

(d)  Del  Giudice  et  al.  [10]  ask  for  clarification  re¬ 
garding  the  time  elapsed  between  resuspension  and 
recording  of  the  spectra.  In  ref.  [9]  we  clearly  speci¬ 
fied  that  “..no  lines  were  ever  detected  in  the  20- 
300  cm-1  region  at  any  time  between  10  min  and 

2  h  30  min  after  resuspension  in  DMM”. 

0.375-9601  /85/S  03.30  ©  Elsevier  Science  Publishers  B.V 
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(e)  All  the  apparent  random  Raman  lines  could  be 
seen  only  in  our  spectra  taken  under  flowing  condi¬ 
tions.  Such  lines  were  narrow  enough  to  be  within 
one  or  two  frequency  steps  (i.e.  with  a  width  of  4—8 
cm-1)  of  the  computer  driven  spectrometer.  This  is 
compatible  with  the  hypothesis  that  either  air  bubbles 
or  microscopic  fluorescent  contaminants  may  have 
caused  these,  by  randomly  passing  through  the  scat¬ 
tering  volume. 

(f)  Under  nonflowing  conditions,  the  spectrum 
of  the  medium  alone  (Davis  Minimal  Medium  1  : 3 

in  H-jO)  was  completely  reproducible,  as  were  the  fea¬ 
tureless  spectra  of  B.  megaterium  suspensions  taken 
at  any  time  after  tesuspension.  It  may  be  useful  to 
point  out  that  all  the  spectra  of  a  single  experimental 
session  were  taken  using  either  the  same  suspension 
for  an  extended  period  of  time  (up  to  45  min)  or  by 
changing  the  suspension  every  10  min  with  another 
sample  taken  from  a  culture  kept  in  a  water  bath  at 
the  selected  temperature. 

Furthermore,  we  have  found  that  Kinoshita  et  al. 
[13]  presented  detailed  calculations  regarding  pos¬ 
sible  artifacts  that  might  have  affected  the  intensity 
ratio  of  Stokes  and  anti-Stokes  lines  in  Raman  spec¬ 
tra  of  Chlorella  pyrenrridosa  [14,15] .  These  possible 
artifacts  include  the  1/A4  dependence  of  the  scatter¬ 
ed  light  intensity,  effective  detection  efficiency,  reso¬ 
nance  effect,  laser-induced  temperature  increase  and 
attenuation  within  the  sample.  These  findings,  along 
with  our  experimental  data,  show  that  several  experi¬ 
mental  artifacts  can  be  introduced  in  the  Raman  spec¬ 
troscopy  of  living  cells  and  that  to  our  knowledge 
no  independent  positive  replication  of  the  data  pre¬ 


sented  in  refs.  [1-3,14,15]  is  available  in  the  litera¬ 
ture. 

In  conclusion  we  do  not  think  that  the  alleged 
presence  of  lines  in  Raman  spectra  of  bateriai  cells 
can  be  claimed  to  be  “..non  questionable...”  [10] 
or  that  such  “...revelations...”  [10]  can  be  taken  as 
a  scientific  proof  of  coherent  phenomena  occurring 
in  living  cells  as  suggested  by  Webb  [3] ,  Frdhiich  [6] 
and  Del  Giudice  et  al.  [7,8,10] . 
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Several  biological  effects  resulting  from  exposure  to  millimeter  waves  have  been  reported.  In  an 
attempt  to  determine  if  millimeter  waves  might  affect  the  conformational  state  of  membranes, 
sonicated  liposomes  have  been  irradiated  with  millimeter  waves  at  41.650  GHz,  stabilized  in 
frequency  to  ±  50  Hz.  Raman  spectra  from  the  lipid  vesicles  were  then  collected  in  the 
conformational  and  in  the  C-H  stretching  region.  No  changes  in  either  the  Raman  peak 
locations  or  relative  intensities  were  detected  upon  millimeter-wave  irradiation,  either  below  or 
above  the  transition  temperature  of  the  phospholipid  (41  *C). 


Several  investigators  have  suggested  that  millimeter 
waves  (mm  waves)  affect  a  variety  of  biological  functions 
including  growth  of  bacterial  cells,12  DNA  and  RNA  syn¬ 
thesis,3  etc.4J  Unfortunately  many  of  these  reports  lack  es¬ 
sential  details  of  the  experimental  procedure,  thereby  severe¬ 
ly  limiting  the  possibility  of  reproducing  the  experiments. 
An  extensive  review  of  the  published  literature  is  available.6 
Many  times  when  duplication  has  been  attempted,  negative 
results  have  been  reported.  In  particular,  no  effects  were  re¬ 
ported  on  cell  morphology7  and  protein  synthesis8  in  mam¬ 
malian  cells.  BHK-21.  and  the  growth  rate  of  E.  coli 9  and 
yeast"’  was  not  altered  by  exposure  to  mm  waves.  In  1968, 
with  the  appearance  of  the  first  experimental  reports.  Froh- 
lich  proposed  that  mm  waves  may  affect  living  systems  by 
inducing  coherent  excitations  in  some  parts  of  cells  like 
DNA.  cell  membranes,  and  enzyme-substrate  complex¬ 


es.4 11  Even  though  no  specifics  were  given,  two  of  the  more 
interesting  features  of  his  model  are  its  frequency  depen¬ 
dence  and  an  energy  threshold  for  the  onset  of  the  effect. 
This  interesting  hypothesis,  however,  has  not  led  to  a  predic¬ 
tion  of  specific  frequencies  or  sites(  s )  of  action  for  it.  Never¬ 
theless,  based  on  the  typical  membrane  thickness  ( 100  A) 
and  on  the  speed  of  sound.  Frohlich  has  suggested  that  oscil¬ 
lations  in  the  10‘°-10"  Hz  range  may  be  excited  in  cell  mem¬ 
branes.4 

Cytoplasmic  membranes  are  an  essential  pan  of  living 
cells.  In  addition  to  acting  as  the  boundary  between  the  in¬ 
ternal  and  the  external  cell  environment,  membranes  also 
perform  essential  functions  such  as  selective  ion  pumping 
and  cell  recognition. 11  As  depicted  by  the  fluid  mosaic  mod¬ 
el,1'  membranes  are  composed  of  lipids  and  proteins  orga¬ 
nized  in  a  bilayer.  70-100  A  thick.  Most  of  the  functional 
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properties  of'  membranes  arise  from  glycoproteins  anchored 
to  the  bilayer  itself  ‘ 4 

To  study  systematically  the  problem  of  identifying  the 
specific  sitei  si  of  possible  interaction  of  mm  waves,  we  used 
Raman  spectroscopy  to  determine  if  mm  waves  can  induce 
conformational  changes  in  sonicated  vesicles  exposed  to 
dl  oiO  GHz.  a  frequency  at  which  relatively  large  effects 
have  been  reported  in  the  past  on  the  growth  rate  of  yeast. ; 

Phospholipids  are  characterized  by  a  temperature  T. 
that  marks  the  phase  transition  between  gel  ( T  <TC )  and 
liquid-crystalline  phase  ( T>  Tc ). 15  In  natural  and  artificial 
membranes  t. like  liposomes),  this  transition  temperature 
and  the  melting  curve  are  known  to  be  affected  by  several 
factors,  including  lipid  moiety1®  and  presence  of  proteins.17 

Dipaimitovl-phosphatidylcholine  (DPPC)  crystals 
( purity  99%  —  )  were  obtained  from  Sigma  and  used  with¬ 
out  further  purification.  For  a  typical  experiment  0.8  mg  of 
DPPC  was  mixed  in  1  ml  of  buffer  (Tris  10  mM.  KC1 0. 1  M, 
pH  7.2 )  and  sonicated  for  30  min.  with  the  tube  containing 
the  suspension  kept  in  an  ice  bath.  The  liposome  suspension 
was  used  after  filtration  through  an  0.8  /xm  Millipore*  filter. 
The  predicted  thickness  of  the  liposomes  is  70  A.1 :  and  the 
average  diameter  is  about  300  A. 18 

Raman  spectroscopy  has  been  widely  used  in  the  past  to 
determine  the  melting  curves  of  phospholipids. 1 4  and  the 
conformational  changes  induced  in  liposomes  by  various 
chemical  and  physical  agents  such  as  sonication.:°  cholester¬ 
ol. :i  and  proteins.'7  For  the  present  study  we  used  Raman 
spectroscopy  to  gather  information  on  possible  effects  of  mm 
waves  on  the  conformation  of  phospholipids  in  liposomes 
above  and  below  Tc . 

The  laser  source  was  a  Spectra  Physics  164/09  argon- 
ion  laser  operating  at  314.3  nm.  A  plasma  line  filter  (Spex) 
was  found  to  be  essential  to  eliminate  nonlasing  plasma  lines 
that  would  otherwise  become  apparent  in  the  Raman  spectra 
arising  from  the  high  scattering  efficiency  of  the  suspension. 
The  power  at  the  sample  location  was  400  mW  Spectra  were 
collected  with  a  Spex  1403  double  monochromator  coupled 
to  a  Spex  third  monochromator.  The  slit  width  was  set  at  400 
pm.  with  a  spectral  bandwidth  of  4  cm  - '  Spectra  were  col¬ 
lected  at  discrete  wave  numbers  with  an  accumulation  time 
of  5  s  and  a  wave-number  step  size  of  3  cm  “ 1  The  Raman 
scattered  light  intensity  was  measured  with  an  RCA  C3 1034 
cooled  GaAs  photomultiplier  connected  to  a  photon 
counter,  Princeton  Applied  Research  1 109  A  custom-made 
microprocessor-based  controller  was  used  to  control  the 
scanning  limits,  integration  time,  and  step  size  of  the  triple 
monochromator  assembly. 
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The  liposome  suspension  was  placed  in  a  borosilicate 
capillary.  0.2  mm  i.d..  that  was  then  placed  in  a  larger  capil¬ 
lary.  0.o5  mm  i  d.,  in  which  temperature  controlled  water 
was  circulated.  The  double  capillary  assembly  was  inverted, 
at  a  45“  angle,  in  a  WR-19  waveguide.  A  diagram  of  the 
waveguide/ capillary  assembly  is  -shown  in  Fig.  I.  A  thermo¬ 
couple  was  placed  inside  the  larger  capillary  to  monitor  ihe 
temperature  of  the  water  jacket.  The  maximum  temperature 
increase  caused  by  the  laser  exposure  in  the  water  jacket  was 
about  2.5  *C. 

The  mm-wave  source  was  a  klystron,  Vanan  302-BT 
The  irradiation  frequency  was  measured  and  phase-locked 
with  an  EIP  578  source-locking  frequency  counter,  with  a 
resulting  short  term  stability  ( 1  h )  of  ~  50  Hz  and  mea¬ 
sured  with  up  to  II  significant  digits.  The  irradiation  fre¬ 
quency  was  4 1.650  GHz  and  was  chosen  on  the  basis  of  being 
a  frequency  at  which  a  relatively  large  effect  was  reported 
previously  on  the  growth  rate  of  veast.7  The  cw  power  in  the 
waveguide  was  30  *  0.5  mW  The  temperature  increase  of 
the  water  jacket  induced  by  the  mm-wave  irradiation  did  not 
exceed  0.2  *C. 

The  power  deposition  pattern  of  the  liposome  suspen¬ 
sion  in  the  capillary  was  estimated  numerically The  calcu¬ 
lations  indicated  that  1 1  %  of  the  incident  power  was  ab¬ 
sorbed  tn  the  sample,  with  a  resulting  average  specific 
absorption  rate  of  about  bOOO  W /kg  The  power  absorbed  by 
the  whole  double  capillary  assembly  was  80%.  This  was  in 
fair  agreement  with  expenmental  data. 

Raman  spectra  were  collected  near  the  conformational 
frequency  shift  region.  700- 1 700  cm'1,  and  C-H  stretching 
region,  2800-3500  cm  ~ 1  To  study  the  possible  effects  of  mm 
waves  on  liposomes  both  in  the  gel  and  m  the  liquid-crystal¬ 
line  phase,  spectra  were  taken  at  25  and  45  ”C.  i.e..  below  and 
above  7;,41  *C,  of  DPPC. 
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Spectra  collection  was  started  within  1  h  after  sonica- 
tion.  The  fluorescence  background  was  partially  quenched 
by  exposing  the  sample  to  the  laser  beam  for  30  mm  before 
starting  the  collection  of  spectra.  Spectra  at  25  *C  were  taken 
first  without  mm  waves,  and  then  with  mm-wave  irradia¬ 
tion.  Millimeter-wave  irradiation  was  interrupted  before 
bringing  the  suspension  to  45  3C. 

The  spectra  taken  with  and  without  irradiation  are 
shown  in  Figs.  2(a)— 2{d)  for  the  conformational  region, 
and  in  Figs.  3(a)— 3(d)  for  the  C-H  stretching  region.  The 
spectra  were  not  smoothed.  In  agreement  with  previous  re¬ 
ports,"  the  temperature-induced  spectral  changes  are  ap¬ 
parent  from  the  changes  in  intensity  of  the  lines  of  the  skel¬ 
etal  optical  modes  at  1064  and  1 128  cm” 1  and  in  the  CH, 
stretching  modes  at  2882  and  2845  cm”1.  On  the  other  hand, 
no  differences  were  apparent  in  either  region,  or  at  either 
temperature,  between  control  and  mm-wave  irradiated  sam¬ 
ples. 

Cell  membranes  have  been  suggested  in  the  Frohlich 
model  as  one  of  the  possible  sites(s)  of  interaction  of  mm 
waves  with  living  systems.4  Liposomes  were  used  in  the  pres¬ 
ent  study  to  provide  a  first  step  in  identifying  the  possible 
site(s)  of  interaction  of  mm-wave  irradiation  using  a  rela¬ 
tively  simple  and  well -characterized  model  of  cell  mem¬ 
branes.  Raman  spectroscopy  was  used  not  only  for  its  capa¬ 
bility  in  detecting  conformational  changes  in  phospholipids, 
but  also  because  it  allowed  this  determination  to  be  per¬ 
formed  simultaneously  to  the  mm-wave  irradiation.  As  Figs. 
2  and  3  show,  there  was  no  evidence  of  effects  induced  by  the 
mm-wave  irradiation.  This  experiment,  while  it  cannot  con¬ 
firm  Frohlich  hypothesis,  cannot  prove  it  wrong  either. 

Even  though  liposomes,  as  prepared  for  the  present 
study,  cannot  mimic  the  functional  complexity  of  actual  cell 
membranes,  they  offer  the  advantage  of  allowing  accurate 


FIG  }  Raman  spectra  ot'  DPPC  sonicated  vesicles  in  the  C-H  stretching 
region  with  and  without  mm  waves,  a  j  25  *C  nomrradiaied;  ( b )  25  ’C  irra- 
Jiated.  a  45  'C  nomrradiaied:  ■  d)  45  *C  irradiated.  The  intensity  scale  is 
relative  to  spectrum  ui  The  other  spectra  have  been  plotted  with  the  fol¬ 
lowing  offset  bi  J0(XX)  counts/ 5  v  c>  bO  (XX)  counts/ 5  s.  d)  ->0  000 
counts/  5  s. 


identification  of  the  system  being  irradiated.  For  future  ex¬ 
periments,  however,  liposomes  could  be  prepared  by  incor¬ 
porating  integral  proteins  that,  among  other  functions,  will 
simulate  the  presence  of  ion  pores,"  a  structure  essential  to 
membrane  function.  The  incorporation  of  proteins,  how¬ 
ever.  greatly  alters  the  conformational  state  of  phospholipids 
in  vesicles  prepared  with  pure  lipids,  while  the  protein  struc¬ 
ture  is  not  altered  when  incorporated  in  a  lipid  bilayer. ' 

Furthermore,  experiments  at  a  lower  frequency,  2.45 
GHz,  have  indicated  that  microwaves  increase  the  perme¬ 
ability  of  erythrocytes  around  Te,  but  not  at  temperature 
different  from  Tc  ,24  Raman  spectroscopy  could  be  used  to 
determine  if  such  changes  in  permeability  can  be  ascribed  to 
conformational  changes  in  the  lipid  bilayer,  when  exposed 
either  to  2.45  GHz  or  to  mm-wave  radiation. 

Many  effects  have  been  ascribed  to  exposure  to  mm 
waves.  Some  investigators  have  even  suggested  the  use  of 
mm  waves  as  a  possible  diagnostic  and  therapeutic  tool  for 
cancer.4  The  present  study  did  not  show  any  effect  on  the 
conformational  state  of  model  membranes  exposed  at  4 1 . 650 
GHz,  either  below  or  above  Tc.  While  this  study  cannot 
exclude  the  possibility  that  some  biological  functions  might 
be  affected  in  nonthermai  ways  be  millimeter- wave  irradia¬ 
tion,  it  adds  to  a  growing  list  of  reports  showing  the  absence 
of  statistically  significant  effect  when  a  variety  of  living  or¬ 
ganisms  are  exposed  to  millimeter-wave  irradiation.*’"10 
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MEASUREMENT  OF  COMPLEX  PERMITTIVITIES  OF  BIOLOGICAL 
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ABSTRACT 

Dielectric  properties  of  biological  materials  have  been  measured  in 
the  frequency  range  of  26.5  to  60  GHz  using  a  modified  Infinite  sample 
method.  In  this  method,  Impedance  transformers  are  used  to  reduce  the 
reflectivity  of  the  infinite  sample  of  biological  material.  The  complex 
reflection  coefficients  were  measured  by  an  automated  measurement  system 
using  three-probe  ref lectoraeter  algorithm.  An  uncertainty  analysis  is 
performed  to  get  an  estimate  of  the  errors  in  the  measured  complex 
permittivities.  The  experimental  results  and  the  Debye  relaxation 
parameters  are  reported  for  ten  percent  NaCi  solution,  human  blood, 
aqueous  solution  of  bovine  serum  albumin  and  cell  suspensions  of 
saccaromyces  cerevisiae. 
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INTRODUCTION 


Knowledge  of  the  dielectric  properties  of  biological  substances  is 
necessary  for  determining  absorption  of  electromagnetic  energy  in 
biological  bodies  and  for  understanding  biophysical  mechanism  of 
interaction  with  electromagnetic  fields.  Several  investigators  have 
reported  sharp  distinct  resonances  in  the  absorption  of  millimeter  waves 
by  various  biological  preparations  [1-3].  Information  on  dielectric 
properties  is  useful  to  study  these  frequency-dependent  effects  and  weak 
nonthermal  interactions  of  biological  systems  with  millimeter-wave 
radiation. 

Previously,  the  dielectric  properties  of  biological  materials  have 
been  measured  below  24  GHz  by  several  investigators  [4-6].  But,  in  the 
millimeter-wave  range,  the  dielectric  data  have  been  sparse  or 
nonexistent  because  of  (1)  limitations  of  conventional  microwave 
techniques  regarding  small  sample  size  associated  with  high  attenuation, 
(2)  requirements  of  precision  measurement  facilities  and  sophisticated 
instrumentation  and  (3)  difficulties  in  eliminating  experimental 
artifacts  (e.g.,  air  bubbles  and  nonuniform  packing  densities  in  cell 
suspensions,  etc.).  The  dielectric  properties  of  water  have  been 
measured  at  34.86  and  70  GHz  by  Grant  [7]  and  Szwamowski  and  Sheppard 
[8],  respectively.  Also,  Edrich  and  Hardy  [9]  have  measured  complex 
permittivities  of  muscle  and  fat  tissues  at  frequencies  from  40  to  90 
GHz.  For  investigations  of  frequency-dependent  absorption,  the 
millimeter-wave  absorption  spectra  of  different  biological  substance, 
(e.g.,  DNA,  RNA,  etc.)  have  been  measured  by  Gandhi,  et  al.  [10]. 
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Because  of  the  ubiquitous  presence  of  water  in  biological  materials 
and  its  high  attenuation  (15-30  dB/mm)  at  millimeter  wavelengths,  a 
layer  of  less  than  0.3  mm  of  biological  substance  absorbs  80  to  90 
percent  of  the  radiation.  So,  the  Roberts  and  Von  Hippel's  shorted-line 
technique  [11]  and  the  microwave  bridge  technique  (12)  would  require 
suhmi llimeter  sample  lengths  and  complicated  sample  holders.  Therefore, 
for  accurate  broadband  measurements,  it  is  undesirable  to  use  these 
conventional  microwave  techniques.  A  review  of  other  microwave 
measurement  techniques  for  dielectric  measurements  of  biological 
materials  is  given  in  [13,  14],  Attempts  have  been  made  to  extend 
quasi-optical  techniques  towards  the  millimeter-wave  region.  At 
wavelengths  above  5  mm,  the  quasi-optical  technique  (namely;  dispersive 
Fourier  transform  spectroscopy)  cannot  be  used  because  of  weak  radiative 
power  from  the  mercury  vapour  lamp  [15].  The  ocher  quasi-optical 
techniaues  (e.g.,  Mach-Zehnder  spectrometer)  are  also  not  particularly 
suitable  for  biological  samples  at  millimeter  wavelengths. 

The  infinite  sample  method  (sample  thickness  >  3-4  skin  depths)  has 
been  used  to  measure  the  dielectric  properties  of  low  dielectric 
constant  and  high  loss  materials  [16,  22].  This  method  obviates  the 
need  for  complicated  sample  holders  and  submillimeter  sample  thicknesses 
because  of  the  use  of  longer  sample  lengths.  But,  for  high  dielectric 
constant  and  high  loss  biological  materials,  this  method  is  very 
inaccurate  because  of  the  need  to  measure  high  voltage  standing  wave 
ratios  (VSWR)  and  small  phase  shifts  (relative  to  short  circuit).  In 
the  modified  infinite  sample  (MIS)  method,  a  quarter  wavelength  (or  odd- 
multiples  of  Quarter  wavelengths)  impedance  transformer  is  used  to 
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reduce  the  reflectivity  of  the  infinite  sample  of  biological  material. 
With  suitable  choice  of  impedance  transformers,  this  method  gives 
accurate  dielectric  measurements  for  biological  materials. 

This  paper  describes  MIS  method  for  dielectric  measurements  of 
biological  materials.  For  measurement  of  complex  reflection 
coefficients,  automated  measurement  systems  have  been  developed  which 
give  reliable  and  reproducible  data  in  the  frequency  range  of  26.5  -  60 
GHz.  The  three-probe  ref lectometer  algorithm  [20]  has  been  used  to 
improve  the  accuracy  of  reflection  coefficient  measurements. 
Uncertainties  in  complex  permittivities  and  errors  in  the  measured 
reflection  coefficient  have  been  discussed  for  the  MIS  method.  Finally, 
the  complex  permittivity  measurements  were  performed  for  various 
biological  materials  including  10  percent  NaC£  solution,  human  blood, 
aqueous  solution  of  bovine  serum  albumin  (BSA)  and  cell  suspensions  of 
saccaromyces  cerevisiae.  The  complex  permittivity  data  were  analyzed 
using  single  relaxation  Debye  dispersion  and  Cole-Cole  model. 

THEORETICAL  BASIS 

Figure  1  shows  the  schematic  diagram  of  the  sample  holder  for  MIS 
method.  It  is  assumed  that  (1)  only  the  dominant  mode  TEjq  propagates 
through  the  filled  sections  of  the  rectangular  waveguide,  (2)  the 
waveguide  section  containing  the  biological  material  is  sufficiently 
long  so  that  it  can  be  considered  as  an  infinite  sample,  (3)  the 
biological  material  is  homogeneous,  and  (4)  the  dielectric  material  for 
the  impedance  transformer  is  homogeneous  and  isotropic. 

The  characteristic  impedance  Zc^  normalized  with  respect  to  the 
characteristic  impedance  of  the  empty  waveguide  ZQ  is  given  by  [22] 
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where  i  »  1  refers  to  the  Impedance  transformer  and  I  »  2  refers  to  the 
biological  material,  is  the  complex  permittivity  of  the  ith 
layer,  Xc  ■  2a  is  the  cutoff  wavelength  of  the  empty  waveguide  and  X  is 
the  free-space  wavelength. 

The  phase  constant  8^  is  given  by  the  relation 

r  2i]/2 

Si  -  (2w/X)  [ej  -  C  X/Xc)^ J  i  -  1,  2  (2) 

The  measured  reflection  coefficient  (r  ■  [r|  eJ  )  at  z  ■  0  is 

related  to  the  normalized  input  impedance  by  the  following  relation 
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using  transmission  line  equations,  the  input  impedance  *in  is  given  by 


Z  -  Z  . 
in  cl 


Zc2  +  -1  Zcl  tan  fBl^ 
Zcl  +  j  Zc2  tan 


(4) 


By  rearranging  various  terms,  we  get 
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From  Eq .  1,  Che  complex  permicclvicy  of  biological  material  is  given  by 


By  substituting  £q.  1-5  in  Eq.  6,  one  can  obtain  an  explicit  expression 

for  the  complex  permittivity  of  the  biological  material  in  terms  of 
» 

T,  (dielectric  constant  of  impedance  transformer)  and  l  (effective 

length  of  the  impedance  transformer).  It  is  assumed  that  the  dielectric 
materials  for  impedance  transformers  are  low-loss  materials  (e.g., 

**  f 

alumina).  So,  c1  is  negligible  compared  with  c^. 

IMPEDANCE  TRANSFORMERS  AND  SAMPLE  HOLDERS 

Figures  2  and  3  show  the  variation  of  |r|  and  0  for  the  infinite 
sample  of  water  with  and  without  alumina  impedance  transformers.  Some 
important  observations  can  be  made  by  inspecting  Figs.  2  and  3.  For  a 
band  of  frequences  close  to  fwin,  the  effects  of  introducing  impedance 
transformers  are  (1)  the  increased  variation  in  |r(  and  0  and,  (2)  the 
reduction  in  |r|  from  the  values  greater  than  0.76  to  less  chan  0.4. 
Because  of  these  reasons,  for  a  given  precision  in  of  water,  it  is 
possible  to  allow  higher  measurement  errors  in  jrj  and  0  by  using 
impedance  transformers.  This  is  shown  in  Table  1  for  several 
frequencies  by  implementing  a  computer  algorithm  using  Eq .  1-5.  Since 
it  is  difficult  to  measure  phase  with  an  accuracy  less  than  0.5°, 
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particularly  at  millimeter  wavelengths,  the  introduction  of  a 
transformer  is  advantageous. 

Another  advantage  of  using  an  impedance  transformer  is  that  it 
obviates  the  need  for  precision  measurements  of  high  VSWRS.  Errors 
Introduced  in  the  measurements  of  high  VSWRs  are  those  due  to  the 
absence  of  strict  square  law  response,  appreciable  rise  and  fall  time 
for  square  wave  source  modulation  and  probe  reflections  due  to  excessive 
probe  penecration  [17-19].  Because  of  the  limitations  of  millimeter 
wave  instrumentation,  many  of  these  errors  cannot  be  reduced 
significantly. 

The  impedance  transformers  may  be  quarter  wavelength  (or  odd- 
multiples  of  auarter  wavelength)  waveguide  sections  of  suitable  solid 
macerlal.  In  order  to  achieve  the  largest  possible  bandwidth  (band  of 
frequencies  where  f  f  (  <  0.4),  it  is  necessary  to  use  a  transformer  of 
length  one-quarter  or  three-quarter  wavelength.  For  higher  odd- 
multiples  of  quarter  wavelengths  there  will  be  rapid  variation 
of  jrj  and  9  which  limits  the  useful  bandwidth.  It  is  required  to 
select  a  material  with  dielectric  constant  between  six  and  ten  because 
of  impedance-matching  considerations.  Alumina  is  chosen  as  the 
impedance  transformer  material  because  of  its  low-loss,  good 
machinabi lity  and  appropriate  dielectric  constant  (t^  ■  9.25  ).  In 

order  to  cover  the  frequency  range  of  26.5  -  60  GHz,  three  alumina 
transformers  of  lengths  0.069,  0.163  and  0.126  cm  were  fabricated  for 
frequency  ranges  28-33  GHz,  42-45  GHz  and  52-57  GHz,  respectively.  It 
should  be  noted  that  one-quarter  wavelength  transformer  is  used  for  28- 
33  GHz  frequency  range,  while  three-quarter  wavelength  transformers  are 


used  for  42— *5  GHz  and  32-57  GHz  frequency  ranges.  Because  of 
difficulties  in  fabrication  of  transformer  thicknesses  less  than  0.05 
cm,  the  three-quarter  wavelength  transformers  were  used  for  higher 
frequency  ranges. 

The  complete  sample  holder  consists  of  the  impedance  transformer 
and  a  section  of  the  waveguide  containing  the  biological  sample.  A 

plastic  (Mylar”,  0.001”  chick)  window  is  used  at  z  ■  L  Interface  (Fig. 
1),  so  as  to  prevent  the  flow  of  biological  sample  into  impedance 

transformer  and  slotted  line.  The  dielectric  discontinuity  due  to 
plastic  window  and  discontinuities  due  to  possible  slight  misalignment 
of  flanges  at  the  interfaces  z  ■  i  and  z  *  0  (Fig.  1)  were  taken  into 
account  by  a  small  change  (<  2  percent)  in  the  length  of  impedance 
transformer.  This  is  Justified  since  small  discontinuities  can  be 
modeled  as  susceptances  which  lengthen  or  shorten  the  transmission 
lines.  The  effective  length  of  the  impedance  transformer  is  measured  by 
using  standard  materials  (e.g.,  water,  saline,  etc.)  as  the  calibrating 
med la. 

THREE-PROBE  REFLECTOMETER 

In  order  to  make  precision  measurement  of  J  r [  and  9  using  an 

automated  measurement  system,  the  three-probe  ref lectometer  formulation 
described  by  Caldecott  [20]  has  been  implemented.  It  uses  standing  wave 
(SW)  voltage  measurements  for  three  equally  spaced  positions  of  the 
slotted  line  probe.  It  is  assumed  that  the  slotted  line  detector  is  in 
the  square  law  region  and  the  frequency  is  known  accurately.  Let  V  be 
the  incident  voltage  and  Vn  is  the  SW  voltage  at  the  nth  probe 
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position.  Then 


vn  -  V  [i  .  |r|  .'<e  •  9n)]  <7) 

pn  ’  lv„|2  •  v2  [>  *  lr|2  ♦  2|rl  =°»  (9  "  «„)] 

where  9  is  the  phase  shift  corresponding  to  the  distance  to  the  load 
n 

and  back.  There  are  three  unknowns  (|r|,  0  and  V  )  in  Eq.  8.  They  are 

calculated  from  the  three  measured  values  of  Pn.  The  complete 

formulation  is  given  in  [20].  The  spacing  between  the  neighboring  probe 

positions  is  close  to  one-eighth  of  a  guide  wavelength  ( Xg/8 )  for 

maximum  accuracy  in  jrj  and  9.  The  j  T J  in  this  algorithm  is  inversely 
2  2 

proportional  to  V  ,  while  6  is  independent  of  V  .  So  the  errors  in  the 

•>  i  | 

calculated  values  of  V"  cause  inaccuracies  in  f  T | •  For  high  reflection 

coefficient  loads  (e.g.,  short  circuit,  infinite  sample  of  water,  etc.) 

2 

the  errors  in  the  calculated  values  of  V  due  to  errors  in  the  measured 
SW  voltages  are  large.  Hence,  high  values  of  ( T  |  cannot  be  measured 
accurately  by  this  algorithm.  Since  reflection  coefficient  measurements 
for  MIS  method  Involve  |  r |  <  0.5,  this  limitation  is  not  important  for 
our  measurements. 

MEASUREMENT  SYSTEM 

Two  computer-controlled  measurement  systems  covering  the  frequency 
ranges  26.5  -  40  GHz  and  40  -  60  GHz  have  been  developed  for  measurement 
of  complex  reflection  coefficient.  The  measurement  system  (40  -  60  GHz) 
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is  schematically  illustrated  in  Fig.  4.  The  IMPATT  oscillators  have  an 
amplitude  and  frequency  accuracy  of  ±  0.06  dB  and  ±  3  MH2  [10].  The 
ferrite  modulator  provides  1  KHz  square  wave  modulation  of  the 
millimeter  wave  signal.  The  travelling  probe  of  the  slotted  line  was 
moved  in  steps  of  0.3  mm  by  the  computer-controlled  stepping  motor,  so 
as  to  take  40  SW  voltage  readings  of  the  standing  wave  pattern.  Because 
of  computer  control  of  the  stepping  motor,  an  accuracy  of  ±  0.01  mm  is 
realized  for  all  probe  positions.  At  each  location  of  the  probe,  the  SW 
voltage  was  measured  by  a  tunable  diode  detector  (sensitivity  *  -  60 
dBm)  via  SWR  meter,  A/D  converter  and  FP-86B  computer.  The  SWR  meter 
was  used  as  a  high  gain  amplifier  for  the  detected  1  KHz  signal.  The  SW 
voltages  are  averaged  over  ten  or  more  repetitive  measurements 
(depending  on  the  variance).  Because  of  averaging,  the  precision  of  SW 
voltages  is  better  than  ±  0.05  dB. 

The  three-probe  ref lectometer  algorithm  [20]  (which  is  implemented 
on  HP-86B  computer)  is  used  to  calculate  more  than  20  sets  (depending  on 
frequency)  of  |  T  (  and  0.  The  accuracy  of  the  complex  reflection 
coefficient  (r  »  f  T  |  e^9ni)  is  improved  by  averaging  the  various  values 
of  T  and  9  .  At  each  frequency,  it  is  necessary  to  ensure  that  the 
diode  detector  is  in  the  square  law  region  because  of  the  changes  in 
power  output  of  the  source  and  the  detector  response  with  frequency. 
This  is  done  by  adjusting  the  power  level  using  the  attenuator.  The 
absolute  phase  9  (degrees)  is  obtained  from  the  measured  phase  9 

IB 

(degrees)  by  estimating  the  length  l  between  the  starting  probe  position 
and  the  load.  The  following  relationship  is  obtained  from  the 

transmission  line  theory. 


L 


.■sAa 


(9) 


0  -  9  +  720  i/Xfr 

w 

The  length  1  is  calculated  from  short  circuit  phase  measurement  by  using 
Eq .  9  for  9  -  180*.  An  accuracy  of  the  order  of  0.0016  Xg  is  realized 
for  length  l  by  frequency  averaging  of  the  length  data  and  repetitive 
short  circuit  measurements.  The  measurements  were  not  made  in  the 
frequency  ranges  35  -  40  GHz  and  46  -  50  GHz  because  of  very  low  probe 
output  from  the  slotted  line  and  amplitude  instability  of  the  IMPATT 
oscillator.  There  were  a  few  frequency  ranges  such  as  33.5  -  34.5  GHz 
and  57.5  -  60  GHz  where  the  measurements  were  not  made  because  of 
frequency  instability  and  low  power  output  of  the  source. 

UNCEKTAINTT  ANALYSIS 

The  total  uncertainty  in  measurement  of  complex  permittivity 
-  z'n  -  je")  usin?  the  -IS  method  mainly  depends  on  (1)  errors  in 
measured  |  ~  |  and  9  ,  and  (2)  the  errors  in  the  measured  dielectric 
constant  and  length  of  the  Impedance  transformer. 


The  total  uncertainties  Ae^  and  Ac"  are  defined  as  [21] 


(11) 
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where  4e'^  and  41  are  the  errors  in  the  dielectric  constant  and  the 

effective  length  1  of  the  impedance  transformer.  A|r|  and  49  are  the 

errors  in  the  magnitude  and  phase  of  Che  reflection  coefficient.  The 
partial  derivatives  appearing  in  Eq.  10  and  11  cannot  be  evaluated 
analytically  for  the  MIS  method.  So,  a  computer  algorithm  using  Eq .  1  - 
6  was  used  to  calculate  and  Ac!]. 

The  primary  sources  of  errors  in  |l*j  and  9  are  the  measurement 
errors  in  SW  voltages  and  inaccuracies  in  the  position  of  the  slotted 
line  probe.  The  other  likely  sources  of  errors  are  the  reflections  from 

the  probe,  amplitude  and  frequency  drift  of  the  IMPATT  oscillators, 

square  law  response  errors,  appreciable  rise  and  fall  time  of  the  square 
wave  source  modulation,  etc.  Based  on  our  measurements  with  short 
circuit  .and  standard  materials  (e.g.,  water),  it  is  estimated  chat  the 
errors  in  J  T  |  and  9  are  less  than  ±  1.5  percent  and  ±  1°, 

respectively.  Table  2  shows  the  uncertainties  of  and  e]j  for  ±  1.5 
percent  and  ±  1°  errors  in  |r|  and  9,  respectively. 

The  dielectric  constant  of  the  impedance  transformer  material 
(Alumina)  was  measured  at  several  frequencies  (30  -  60  GHz)  using  the 
shorted  line  technique  described  by  Roberts  and  Von  Hippel  [11].  The 
measured  data  give  *  9.25  *  0.05.  This  value  is  in  fair  agreement 
with  those  in  the  literature  [23,  15]  considering  the  differences  in 

manufacturers,  types  of  alumina  and  frequency  ranges. 

The  effective  lengths  of  impedance  transformers  are  estimated  by 

—A 

using  water  having  a  conductivity  of  200  uS/m  (a  •  2  x  10  S/m)  as  a 
standard  material.  Debye  parameters  for  water  at  23®C  are  assumed  as 

c  -  78.56,  t  ■  4.3  and  t  ■  8.58  ps  [7],  where  £  ,  t  and  T  are  the  low 
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frequency  permittivity,  high  frequency  permittivity  and  relaxation  time, 
respectively.  The  physical  lengths  i  are  varied  (within  ±  2  percent), 
so  as  to  obtain  minimum  root -mean-square  error  (RMSE).  The  RUSE  is 
defined  as 


1  r  |  ,,  ..  df.  >21/2 

2N  |£2  (fi,£)  ~  E  ^fi^| 

a.  m 


where  (f^*  is  the  measured  complex  permittivity  calculated  using 
|r|  ,  8,  f^,  ej  ,  and  1.  is  the  ith  measurement  frequency.  ls 
Che  complex  permittivity  calculated  using  single  relaxation  Debye 
dispersion  model.  N  is  the  number  of  measurement  frequencies.  The 
c^fj)  is  defined  as 


where  c  is  the  free  space  permittivity.  Figure  5  shows  t  (f  ,£  ) 

o  Z  1  mi  n 

and  e,{f.)  for  various  frequencies,  where  £  is  the  estimated 
*  l  min 

effective  length  corresponding  to  minimum  RMSE.  The  calculated  RMSE  of 

0.78  for  water  gives  an  estimated  accuracy  of  *  2  urn  for  the  effective 

lengths  of  impedance  transformers.  Table  3  shows  the  physical  lengths 

and  estimated  effective  lengths  (t  )  for  impedance  transformers. 

min 

Table  4  shows  Che  uncertainties  in  and  e“  because  of  ±  0.54 
percent  and  ±  2  urn  errors  in  the  dielectric  constant  and  the  effective 
lengths  of  impedance  transformers.  The  total  uncertainties  in 
and  e“  calculated  using  Eq .  10  and  11  and  Table  2  and  4  are  less  than  ± 
9.5  percent  and  ±  4.5  percent  for  the  frequency  range  of  26.5  to  60 
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GHz.  In  general,  the  errors  in  the  real  and  imaginary  part  of  complex 
permittivities  are  less  than  these  estimated  uncertainties. 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 

In  the  absence  of  any  f reauency-dependent  resonant  effects,  the 
dielectric  properties  of  biological  materials  are  characterized  by  a 
single  relaxation  Debye  dispersion  or  the  Cole-Cole  equation  because  of 
the  relaxation  of  water  molecules.  At  millimeter  wavelengths,  the 
contributions  due  to  other  known  relaxations  (e.g.,  relaxation  of  the 
bound  water,  protein  molecules,  etc.)  are  negligible.  The  Cole-Cole 
equation  is  given  by 


while  the  single  relaxation  Debye  dispersion  is  given  by  Eq.  13.  The 
Cole-Cole  parameter  a  indicates  a  spread  of  relaxation  time  centered 
about  T.  The  complex  permittivity  data  are  fitted  to  Cole-Cole 

equation  and  single  relaxation  Debye  dispersion  by  the  complex  least- 
square  formulation  similar  to  the  one  described  by  Sheppard,  et  al. 
f  24  ] .  For  all  biological  material  used  in  our  research,  it  is  found 
that  the  RUSE  is  not  significantly  different  for  single  Debye  relaxation 
model  (a  ■  0)  and  Cole-Cole  model.  So,  the  single  relaxation  Debye 
dispersion  model  is  adequate  for  the  analysis  of  our  complex 
permittivity  data.  All  the  measurements  were  made  at  the  room 
temperature.  But,  the  maximum  change  in  temperature  (from  the  mean 
temperature)  during  the  experiments  was  limited  to  ±  1®C. 


Figures  6-10  represent  experimentally-measured  and  calculated  (from 
Debye  dispersion)  complex  permittivities  for  10  percent  NaCZ  solution, 
human  blood,  aqueous  solution  of  BSA  and  cell  suspensions  of 
saccarooyces  cerevisiae  for  two  different  concentrations.  Table  5  shows 
the  Debye  parameters  along  with  their  standard  deviations,  mean 
temperatures  and  RflSE  for  various  biological  materials.  For  10 
percent  NaCI  solution,  the  ionic  conductivity  is  calculated  from  Eq.  12 
of  [25].  The  ionic  conductivity  for  whole  human  blood  is  taken  as  1.2 
S/m  from  Cook's  results  [6].  It  is  assumed  that  these  static  values  of 
conductivity  are  valid  at  high  frequencies  because  of  the  negligible 
variation  of  conductivity  with  frequency  for  concentrated  electrolytic 
solutions  [M.  For  all  ocher  biological  materials,  the  ionic 
conductivity  is  taken  as  zero  because  the  measured  static  conductivities 
(using  conductivity  meter)  were  less  than  1  mS/m. 

Figure  6  shows  the  measured  and  calculated  complex  permittivities 
for  10  percent  NaCZ  solution.  Stogryn  [25]  has  given  equations  for 
calculations  of  Debye  parameters  using  previously-published  data. 
The  Cm  value  has  been  assumed  to  be  4.9.  For  10  percent  NaCZ  solution 
at  23*C  c  ^ ,  o  and  T  calculated  from  Stogryn's  equations  are  52.2,  13.4 
S/m  and  7.5  ps,  respectively. 

Considering  the  differences  in  the  value  used  by  Stogryn  and  the 
value  estimated  from  our  data  (See  Table  5),  the  agreement  in  and  T 

is  reasonable. 

The  measured  complex  permittivities  of  whole  human  blood  are  given 
in  Fig.  7  along  with  the  microwave  frequency  data  obtained  from  Cook 
[6].  The  Cook's  data  at  23*C  were  obtained  by  linear  interpolation  of 
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che  data  at  15,  25,  and  35°C.  The  single  relaxation  Debye  dispersion 
equation  was  fitted  to  our  millimeter-wave  data  and  Cook's  microwave 
frequency  data.  These  Debye  parameters  are  given  in  Table  5.  From  [6], 
the  Debye  parameters  for  human  blood  at  25*C  are 

e  *  58.0,  c  ■  4.5,  t  -  9.0  ps  and  c  -  1.2  S/m. 

s  • 

Allowing  for  lower  temperatures  and  possible  differences  in  blood 
samples  regarding  cell  concentration,  these  Debye  parameters  compare 
favorably  with  those  given  in  Table  5.  The  blood  sample  was  taken  from 
a  healthy  male  and  the  red  blood  cell  concentration  was  between  4.8  to 
5.0  x  106  per  mm^.  2.0  mg  of  heparin  was  added  to  the  10  mi  sample  so 
as  to  prevent  coagulation.  In  addition  to  stirring  the  sample,  the 
blood  in  the  sample  holder  was  replaced  by  the  fresh  blood  (from  the 
blood  sample)  auite  frequently  during  the  experiment.  This  reduces  the 
effect  of  nonuniform  cell  packing  density  and  heterogeneity.  The  same 
care  was  taken  while  using  saccaromvces  cerevisiae  cell  suspensions. 

Figure  8  shows  the  measured  and  calculated  complex  permittivities 
for  BSA  solution.  The  solution  was  prepared  at  a  concentration  of  200 
mg/m l  (20  percent)  by  dissolving  BSA  in  deionized  water.  The  BSA  used 
was  supplied  by  CALBIOCHEM.  The  pH  of  the  solution  is  measured  to  be 
fc.9.  At  microwave  frequencies,  the  dielectric  properties  of  7  percent 
BSA  solution  have  been  reported  by  Buchanan,  et  al.  [5].  But,  no 
published  data  are  available  for  20  percent  BSA  solution  at  microwave 
and  millimeter-wave  frequencies.  Because  the  protein  content  of  many 
biological  tissues  is  close  to  20  percent,  it  is  of  interest  to  know  the 


48  - 


dielectric  behavior  of  20  percent  BSA  solution.  The  Debye  parameters 
obtained  from  the  measured  complex  permittivities  are  given  in  Table  5. 

The  protein  solutions  at  millimeter  wave  frequencies  behave  as  a 
nonpolar  solute  (hydrated  protein  molecule)  in  a  polar  solvent 
(water).  So,  it  is  possible  to  calculate  volume  fraction  (V  )  of  the 
hydrated  protein  molecule  by  using  an  appropriate  mixture  equation.  In 
the  literature  the  Maxwell-Fricke  mixture  equation  has  been  employed 
quite  frequently  for  such  purposes  [26,  27],  So,  we  have  used  this 
mixture  equation  for  calculating  the  volume  fraction  of  the  hydrated 
protein  molecule.  The  Maxwell-Fricke  equation  is  given  by 


where  t  ,  e  and  e  are  the  relative  permittivities  of  the  mixture, 
m  w  p 

water  and  hydrated  protein  molecule,  respectively.  This  equation  is 

strictly  valid  for  suspensions  of  spherical  particles.  For  non- 

spherical  particles  the  factor  2  is  replaced  by  a  function  of  the  shape 

of  the  particle-.  In  case  of  BSA,  the  shape  of  the  molecule  is  a  prolate 

ellipsoid  with  an  axial  ratio  of  3  [30].  The  change  in  factor  of  2  is 

small  for  BSA  molecule  because  of  the  small  variation  of  this  factor  for 

prolate  ellipsoids  of  axial  ratio  less  than  5  [31].  From  Table  5, 

c  is  58.5  while  the  c  is  assumed  to  be  5  [26].  e  is  taken  as  78.56 
m  p  w 

from  [7],  So,  Vp  -  0.206.  The  bound  water  can  be  calculated  from  Vp  by 
the  following  relation  [28], 

V  -  C  (V  ♦  5)  (16) 

P 
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where  C  Is  the  concentration  of  protein  solution  in  gms  per  100  mi 
and  V  is  Che  partial  specific  volume  which  is  0.73  for  BSA  molecules 
[30].  From  Eq .  16  the  bound  water  is  calculated  as  0.302  g/g  of 
proteih.  This  is  in  good  agreement  with  the  value  of  0.3  obtained  by 
Buchanan,  et  al.  and  0.42  obtained  from  Viscosity  measurements  [29]. 

Figures  9  and  10  give  the  measured  and  calculated  complex 
permittivities  for  10  and  10  cells/  mi  saccaromyces  cerevisiae  cell 
suspensions.  The  cells  used  in  this  study  were  grown  on  sabouraud 
glucose  agar  (Neopetone**  10  g/i,  glucose  40  g/i,  Baeto  agar**  15  g/i)  and 

stored  at  4®C.  From  these  cultures,  organisms  were  inoculated  in 

R 

sabouraud  glucose  broth  (Neopetrone  10  g/i,  glucose  20  g/i,  ph  5.6), 
incubated  at  32°C  on  an  orbital  shaker  (40  -  50  rpm)  and  centrifuged  in 
the  table  top  centrifuge.  No  published  data  are  available  on  the 
dielectric  properties  of  saccaromyces  cerevisiae  cells  at  microwave  and 
mi llineter-wave  frequencies.  The  volume  fraction  of  the  cellular 


material  and  associated  bound  water  can  be  calculated  from  Eq.  15  for 
10'  and  10^  cells/mi  suspensions.  The  permittivity  (Ep)  Is  assumed  to 
be  5.  This  assumption  is  justified  because  of  the  nonpolar  nature  of 


the  cellular  material  and  associated  bound  water  at  millimeter-wave 

7  Q 

frequencies.  From  Table  5,  t  values  for  10  and  10  cells/mi  are  78.3 

n 

and  61.04.  So,  the  volume  fractions  for  10^  and  10^  cells/mi.  samples 


are  calculated  (from  Eq.  15)  as  0.003  and  0.178,  respectively. 


CONCLUSIONS 

A  new  method  has  been  developed  for  measuring  complex 
permittivities  of  biological  materials  at  millimeter  wavelengths.  This 


method  employs  an  Impedance  transformer  to  reduce  the  reflectivity  of 
che  infinite  sample  of  the  biological  material.  Recause  of  automated 
measurement  system,  good  accuracy  and  reproducibility  were  observed  in 
the  measured  complex  reflection  coefficient  data.  The  experimental 
artifacts  (e.g.,  air  bubbles,  nonuniform  packing  density,  etc.)  were 
minimized  by  the  careful  design  of  experimental  procedures.  From 
uncertainty  analysis,  it  is  estimated  that  the  maximum  errors  in  real 
and  imaginary  part  of  complex  permittivities  are  less  than  9.5  percent 
and  4.5  percent,  respectively. 

The  results  for  10  percent  NaCi  solution  and  whole  human  blood  were 
compared  with  the  previously  published  data.  Good  agreement  was 
observed  regarding  Debye  relaxation  parameters  between  our  data  and  the 
published  data.  The  amount  of  bound  water  was  calculated  for  the 
aqueous  solution  of  BSA.  It  compares  favorahly  with  the  value  reported 
hv  other  investigators.  Finally,  the  dielectric  measurements  were 

1  Q 

performed  on  10'  cells/ml  and  L0  cells/mZ  saccaromyces  cerevisiae 
suspensions.  By  using  Maxvell-Fricke  mixture  eauation,  the  volume 
fractions  were  calculated  for  these  two  concentrations.  The  complex 
permittivity  data  from  all  the  biological  samples  can  be  explained  on 
the  basis  of  Debye  relaxation  of  the  water  molecule  and  no  frequency- 
dependent  effects  were  observed. 
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Table  3.  Physical  and  effective  lengths 
of  Impedance  transformers. 
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Table  5.  Debye  parameters  of  biological  materials 
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Temperature 
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i  c 

Material  1 
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10  percent 
NaCi 

solution 


Human  blood 


BSA  solution 
200  mg/mi  . 


S.  cerevisiae 
107  cells/tni 


S.  cerevisiae 
10*  cells/mi 


55.9*0.9 

7.1*0.13 

4.35*0.25 

58.0±1.6 
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' 
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58.5*1.4 

7.80*0.19  | 
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impedance  transformer  for  an  infinite 
sample  of  water. 


Block  diagram  of  automated  measurement  system 


Fig.  7.  Measured  complex  permittivities  (  •  real,  *  imaginary) 
and  complex  permittivities  (solid  line)  calculated 
from  Debye  dispersion  for  human  blood.  Results 
below  25  GHz  are  taken  from  (6]. 


Fig.  9.  Measured  complex  permittivities  (  •  real,  *  imaginary) 
and  complex  permittivities  (solid  line)  calculated 
from  Debye  dispersion  for  107  cells/  mi 
saccaromyces  cerevlslae  sample. 
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Fig.  10.  Measured  complex  permittivities  (o  real,  *  imaginary)  and 
complex  permittivities  (solid  line)  calculated  from  Debye 
dispersion  for  109  cells/m£  saccaromyces  cerevisiae  sample 
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ABSTRACT 

A  computer-controlled  measurement  system  has  been  developed  for 
complex  permittivity  measurements  of  human  skin  in  vivo  at  millimeter 
wavelengths.  An  area  of  the  human  body  is  pressed  against  the  waveguide 
flange  for  measurement  of  complex  reflection  coefficient.  A  quarter 
wavelength  impedance  transformer  is  used  to  reduce  the  reflectivity  of 
the  human  skin.  The  complex  permittivity  of  the  human  skin  in  vivo  was 
calculated  from  the  measured  complex  reflection  coefficient  by  using  a 
variational  formulation  involving  multimode  analysis.  A  single 
relaxation  Debye  dispersion  is  fitted  to  our  millimeter-wave  data  and 
the  microwave  freauency  data  reported  by  other  investigators.  The 
complex  permittivity  values  for  human  skin  in  vivo  were  much  smaller 
chan  those  reported  for  excised  human  skin  tissue  used  in  the  past.  The 
difference  in  complex  permittivity  values  between  excised  skin  tissue 
and  skin  tissue  in  vivo  may  be  due  to  higher  solid  content  of  the  skin 
in  vivo. 


*This  work  was  supported  by  the  US  Air  Force  School  of  Aerospace  Medi¬ 
cine,  Brooks  Air  Force  Base,  Texas,  under  contract  F  336 1 5-84-K-063 1 . 
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INTRODUCTION 


With  the  increasing  use  of  millimeter-wave  radiation  for 
communications  and  radar,  it  is  necessary  to  know  the  dielectric 
properties  of  various  human  tissues  for  evaluation  of  potential  hazards 
to  humans.  At  millimeter-wave  frequencies,  the  absorption  of 

electromagnetic  radiation  is  mostly  restricted  to  the  skin  because  of 
the  subnillimeter  depths  of  penetration  [1].  So,  the  knowledge  of 
dielectric  properties  of  human  skin  in  vivo  is  of  prime  importance  for 
quantifying  hazardous  effects  of  electromagnetic  waves  in  the 
millimeter-wave  range.  Also,  for  the  study  of  millimeter-wave 
thermography,  it  is  useful  to  know  the  dielectric  properties  of  live 
human  skin.  In  the  past,  the  dielectric  properties  of  excised  human 
skin  tissue  have  been  measured  below  30  GHz  by  Cook  [2],  England  [3,4] 
and  Schwan  [5,  6].  Since  the  dielectric  properties  of  biological 
tissues  above  100  MHz  are  determined  by  their  water,  salt  and  protein 
concent  [7],  the  use  of  dielectric  properties  of  excised  skin  samples  to 
characterize  live  skin  is  justified  only  if  there  is  no  difference 
between  excised  and  live  tissue  materials  regarding  their  water,  salt 
and  protein  concent.  Hence,  the  variation  of  blood  content  between  live 
and  excised  skin  tissue  can  cause  significant  changes  in  the  complex 
permittivities.  Because  of  these  reasons,  it  is  desirable  to  measure 
skin  dielectric  properties  under  conditions  in  vivo. 

Recently,  Hey-Shipton,  et  al.  [8]  have  reported  complex 
permittivities  of  human  skin  in  vivo  over  Che  frequency  range  of  8  to  18 
GHz.  But  their  complex  permittivity  values  are  inaccurate  below  11  GHz 
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because  of  reflections  from  the  back,  of  the  skin  at  lower  frecuencies. 


We  have  previously  presented  complex  permittivity  measurements  of  human 
skin  in  vivo  over  the  frequency  range  of  28  to  33  GHz  [9].  In  this 

paper,  we  are  reporting  complex  permittivity  data  for  human  skin  in  vivo 

(palm  region)  between  28-57  GHz.  The  measurements  were  made  on  a  newly- 
developed  computer-controlled  measurement  system  which  measures  the 
magnitude  and  phase  of  the  reflection  coefficient.  The  complex 
permittivities  of  the  skin  are  obtained  from  measured  reflection 
coefficients  by  using  an  algorithm  which  is  an  extension  of  multimode 
analysis  given  by  Galejs  [10]. 

The  Debye  relaxation  model  is  fitted  in  the  least-square  sense  on 
our  millimeter-wave  data  and  microwave  frequency  data  from  Hey-Shipton, 
et  al.  [8].  The  emissivity  of  the  human  body  and  the  attenuation 
coefficient  of  the  human  skin  are  calculated  over  the  frequency  range  of 

10  to  60  GHz  using  the  complex  permittivity  values  obtained  from  the 

Debye  relaxation  model.  The  complex  permittivity  values  of  human  skin 
in  vivo  given  by  the  Debye  relaxation  model  fitted  to  our  data  and  Hey- 
Shipton  's  microwave  data  are  considerably  lower  than  those  given  by 
Johnson  and  Guy  [11]  and  Cook  [2].  For  example,  at  2.45  GHz  the  complex 
permittivity  value  obtained  from  Johnson  and  Guy  [11]  is  (47  -  J16.2) 
for  excised  skin  tissue  while  the  complex  permittivity  value  given  by 
our  Debye  relaxation  model  for  human  skin  in  vivo  is  (25.0  -  j 13.6). 
The  difference  in  the  complex  permittivity  values  between  excised  skin 
tissue  and  human  skin  in  vivo  may  be  due  to  higher  solid  content 
(proteins  and  lipids)  [7],  approximately  53.12  percent  for  skin  tissue 
in  vivo  compared  to  36.56  percent  for  the  excised  skin  tissue. 
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MEASUREMENT  SYSTEM 


A  solid-stace,  computer-controlled  measurement  system  operating  at 
frequencies  from  26.5  to  60  GHz  has  been  developed  for  measurement  of 
magnitude  and  phase  of  the  reflection  coefficient.  The  measurement 
system  is  schematically  illustrated  in  Fig.  1.  The  key  components  of 
the  measurement  system  are  the  precision  slotted  line,  the  diode 
detector  for  standing  wave  measurements,  frequency  counter  and  HP  86B 

computer  for  overall  automation.  The  traveling  probe  of  the  slotted 
line  is  moved  in  steps  of  0.3  mm  hy  the  computer-controlled  stepping 
motor.  A  large  dial  indicator  calibrated  in  0.01  mm  increments  gives 
the  position  of  the  probe.  At  each  location  of  the  probe,  the  standing 

wave  voltage  is  measured  by  the  diode  detector  (sensitivity  -  -  60  dBm) 

via  SWR  meter,  A/D  converter  and  HP-86B  computer.  By  averaging  10  or 
more  repetitive  measurements  of  standing  wave  voltages,  the  extreme 
variation  was  limited  to  ±  0.05  dB.  The  standing  wave  voltages  were 
measured  at  40  equally  spaced  positions  of  the  slotted  line  probe.  The 
phase  and  the  magnitude  of  the  reflection  coefficients  are  obtained  by 
implementing  a  three-probe  ref lectometer  algorithm  given  by  Caldecoct 

[12J.  For  a  given  setting,  between  20  to  35  independent  values  of  the 
phase  and  magnitude  of  Che  reflection  coefficient  are  obtained  due  to 
different  combinations  of  the  standing  wave  voltages.  The  accuracy  of 
Che  reflection  coefficient  measurements  is  improved  by  a  factor  of  5  -  6 
upon  averaging  these  reflection  coefficient  values.  The  absolute 
phase  ♦  (degrees)  of  the  reflection  coefficient  is  obtained  from  the 
measured  phase  $  (degrees)  of  the  human  skin  and  the  estimated  length 
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i  between  the  starting  probe  position  and  the  load  by  using  the 
following  equation! 

♦  •  ♦  +  7201/Xg  (1) 

m 

The  length  l  is  estimated  from  short  circuit  phase  measurements 
while  Xg,  the  guide  wavelength  is  calculated  from  freauency 
measurements. 

SAMPLE  HOLDER  FOR  IN  VIVO  MEASUREMENTS 

Figure  2  gives  an  arrangement  for  measurements  of  the  human  skin 
in  vivo.  An  area  of  the  human  body  is  pressed  against  the  waveguide 
flange  for  measurement  of  reflection  coefficient.  A  quarter  wavelength 
impedance  transformer  of  Teflon"  is  used  to  reduce  the  reflectivity  of 
the  human  skin.  The  palm  region  of  the  human  body  is  chosen  for 
measurements  because  the  skin  is  thick  in  this  area  C  3  mm)  and  has  a 
relatively  smooch  surface. 

CALCULATION  OF  COMPLEX  PERMITTIVITIES 

The  flow  chart  for  calculation  of  complex  permittivities  of  human 

skin  in  vivo  from  Che  measured  complex  reflection  coefficient  is  given 

in  Fig.  3.  The  key  parts  of  this  flow  chart  are  the  calculation  of 

input  admittance,  admittance  transformation  for  Teflon"  impedance 

* 

transformer  and  calculation  of  new  e  guess  by  zero  finding  technique. 

It  is  assumed  chat  human  skin  is  nonmagnetic  and  its  inhomogene¬ 
ities  are  small  compared  with  the  wavelength  of  the  electromagnetic 
waves.  In  order  to  model  live  human  skin  as  semi-infinite  half  space, 
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it  is  necessary  to  assume  that  skin  is  lossy  enough  that  the  reflections 

from  the  back  of  the  skin  layer  can  be  neglected.  In  practice,  these 

assumptions  are  reasonable  above  11  GHz.  The  input  admittance  due  to 

human  skin  in  vivo  can  be  calculated  by  solving  the  radiation  problem  of 

a  rectangular  waveguide  with  the  semi-infinite  flange  terminated  by  a 

semi-infinite  half  space  of  lossy  dielectrics.  We  have  implemented  a 

formulation  given  by  Galejs  [10]  for  stratified  plasma  for  calculation 

of  input  admittance.  This  formulation  gives  the  input  admittance  of  a 

smaller  waveguide  radiating  into  a  bigger  waveguide  containing  the 

biological  material.  For  lossy  dielectrics  (e.g.,  human  skin  in  vivo) 

it  has  been  found  chat  the  semi-infinite  half-space  is  simulated  by  the 

bigger  waveguide  if  the  ratio  of  bigger  to  smaller  waveguide  is  greater 

chan  30.  The  number  of  propagating  and  nonpropagating  TE  and  TM  modes 

in  bigger  and  smaller  waveguides  is  increased  until  the  convergence  is 

obtained  in  input  admittance  values.  It  has  been  found  that  the 

convereed  values  of  input  admittance  are  obtained  if  all  possible  values 

of  m  and  n  (for  TE__  and  Til  modes)  are  taken  for  0  <  n  <  250 

nn  inn 

and  C  <  n  <  250, 

The  admittance  transformation  due  to  Teflon"  impedance  transformer 

is  calculated  using  transmission  line  eouations.  The  implementation  of 

* 

zero  finding  technique  involved  prediction  of  a  new  guess  for  e  which 
corresponds  to  a  zero  of  error  function  s.  The  Muller  Method  with 
Deflation  [13]  is  used  for  calculation  of  zeros  of  error  function  s. 

In  order  to  test  our  formulation  for  calculation  of  input 
admittance,  we  have  compared  our  results  for  the  case  shown  in  Fig.  4  to 
those  obtained  from  the  Waveguide  Handbook  [14],  The  input  admittance 


(0.5  +  j 0 . 167)  calculated  using  our  program  compares  favorably  with  the 
input  admittance  (0.5  +  jO.17)  given  by  the  waveguide  handbook. 

A  flow  chart  similar  to  Fig.  3  has  been  implemented  by  Decreton 
and  Gardiol  [15]  for  nondestructive  testing  of  dielectric  materials. 
From  Fig.  4  of  [15],  the  calculated  complex  permittivity  value  is  (7.2  - 
j5.9)  for  the  following  case: 

VSWR  -  4.0,  Phase  Shift  -  -4° 

Frequency  *  10  GHz,  WR  -  90  waveguide 

The  schematic  diagram  for  this  case  is  similar  to  Fig.  2.  In  the 
Decreton  and  Gardiol  [15]  method,  however,  no  impedance  transformer  was 
used. 

For  this  case,  our  method  gives  the  complex  permittivity  value  of 
(7.6  -  j5.6).  There  is  a  close  agreement  between  these  two  complex 
permittivity  values.  So,  it  provides  a  good  test  case  for  the  accuracy 
of  our  in  vivo  computer  algorithm.  The  standard  materials  (materials  of 
known  dielectric  properties)  with  complex  permittivity  in  the  range  of 
human  skin  in  vivo  are  not  available  for  frequency  band  26.5  -  60  GHz. 
So,  the  in  vivo  computer  algorithm  has  not  been  tested  by  measuring 
complex  permittivities  of  standard  materials. 

EXPERIMENTAL  RESULTS 

In  order  to  calculate  complex  permittivities  using  the  flow  chart 
given  in  Fig.  3,  it  is  necessary  to  know  the  dielectric  constant  and  the 
length  of  the  impedance  transformers.  At  millimeter  wavelengths,  the 


dielectric  constant  of  Teflon™  has  been  found  to  be  2.05  with  an 
accuracy  of  ±  1  percent  [16,  17].  The  lengths  of  the  impedance 
transformers  were  measured  with  an  estimated  accuracy  of  ±  5  urn  from 
electrical  length  measurements.  The  shorted  line  technique  described  by 
Roberts  and  Von  Hippel  [18]  was  used  for  electrical  length  measurements. 

Using  the  measurement  system  illustrated  in  Fig.  1,  the  magnitude 

and  the  phase  of  the  reflection  coefficients  were  measured  for  the 

arrangement  shown  in  Fig.  2.  The  frequency  ranges  28  -  33  GHz  and  42  - 

57  GHz  are  covered  by  using  two  quarter  wavelength  transformers.  The 

measurements  were  not  performed  in  the  freauency  ranges  35  -  40  GHz  and 

46  -  50  GHz  because  of  the  low  probe  output  from  the  slotted  line  and 

power  instability  of  the  IMPATT  oscillator,  respectively.  Based  on  our 

measurements  with  standard  materials  (e.g.,  water)  used  for  modified 

infinite  sample  method  [19],  it  is  estimated  that  the  errors  in  the 

magnitude  and  the  phase  of  the  reflection  coefficients  are  less  than  ± 

1.5  percent  and  ±  1°,  respectively.  The  uncertainties  in  the  calculated 

* 

complex  permittivities  (c  =  e'  -je“)  are  mainly  due  to:  1.  The  errors 
in  measured  reflection  coefficient  (T  *  |  F j  e^)f  and  2.  The  errors  in 
the  dielectric  constant  (e  )  and  the  length  (4)  of  the  impedance 
transformer.  The  procedure  for  calculation  of  uncertainties  in  the 
complex  permittivities  is  similar  to  the  one  described  in  [19,  20].  The 
maximum  uncertainty  in  e'  and  e"  of  the  human  skin  in  vivo  due  to 
measurement  errors  in  |r|,  9,  and  4  is  less  than  14  percent  and  6 
percent,  respectively.  The  choice  of  the  palm  region  of  the  human  body 
for  complex  permictivitv  measurements  minimizes  errors  due  to  improper 
skin  contact. 
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Figure  5  shows  the  measured  complex  permittivity  (•  for  real  part 
and  *  for  imaginary  part)  of  the  first  author's  skin  at  millimeter 
wavelengths.  The  microwave  frequency  data  from  Hey-Shipton,  et  al.  [8] 
are  also  included  in  Fig.  5.  The  complex  permittivity  data  were  fitted 
to  a  single  relaxation  Debye  dispersion;  i.e., 


where  e  and  e  are  the  low-  and  high-frequency  limits  of  the 

s  * 

permittivity,  t  is  the  relaxation  time,  f  is  the  frequency  and  a  is  the 
ionic  conductivity.  For  human  skin  in  vivo,  o  is  assumed  to  be  1.4  S/mr 
from  Cook's  [2]  results  for  excised  human  skin  tissue.  The  Debye 
parameters  (cs>  cm,  and  t)  were  found  using  the  least-sauare 
minimization  given  by  Sheppard,  et  al.  [21].  The  solid  line  of  Fig.  5 
is  calculated  from  Ea .  2  using  the  least-sauare  fitted  Debye  parameters. 

The  important  electrical  parameters  which  determine  the 
interaction  of  electromagnetic  radiation  with  human  body,  the 
attenuation  coefficient  and  the  emissivity,  are  shown  in  Figs.  6  and  7, 
respectively.  Table  1  gives  the  Debye  parameters  for  human  skin  in  vivo 
from  our  measurement  data  (including  8-18  GHz  data  from  Hey-Shipton, 
et  al.  [8])  and  the  Debye  parameters  for  excised  skin  tissue  from  Cook 
[2].  In  Figs.  6  and  7,  the  dashed  line  represents  the  calculated 
parameters  from  Debye  dispersion  for  human  skin  in  vivo  while  the  solid 
line  represents  the  calculated  parameters  from  the  Debye  dispersion  for 
excised  skin  tissue.  Figure  8  shows  the  complex  permittivity  values 
calculated  for  human  skin  in  vivo  and  excised  skin  tissue  from  their 


Debye  dispersions.  The  root-mean-square  error  (R.MSE)  for  Che 
experimental  data  of  Fi g.  5  has  been  found  to  be  0.62.  The  R MSE  is 


where  is  the  measured  complex  permittivity  at  ith  frequency,  is 
the  calculated  complex  permittivity  from  single  relaxation  Debye 
dispersion  at  the  ith  frequency  and  N  is  the  number  of  data  points. 

DISCUSSION  AND  CONCLUSIONS 

From  Fig.  8,  it  is  observed  that  e'  and  e"  values  for  human  skin 
in  vivo  are  smaller  than  those  obtained  from  excised  human  skin 
tissue.  These  differences  in  permittivity  values  result  in  significant 
changes  in  emissivity  and  attenuation  coefficient  at  millimeter  wave¬ 
lengths.  As  shown  in  Fig.  6,  the  attenuation  coefficients  (a(dB/mm)) 
for  human  skin  in  vivo  are  lower  chan  those  obtained  for  the  excised 

skin  tissue.  So,  the  skin  depth  (-  8.686  /a)  for  live  skin  is  higher 
Chan  that  for  the  excised  skin  tissue.  The  emissivity  (Z)  represents 

the  percentage  of  power  absorbed  by  the  skin.  From  the  results 
presented  in  Fig.  7,  it  is  observed  chat  the  power  coupling  coefficient 
to  Che  human  body  is  somewhat  higher  than  previously  thought  from  Cook 
[2]  and  England  [4]  data  for  excised  skin  tissue. 

In  Table  1,  Che  Dehye  parameters  for  excised  skin  tissue  given  by 
Cook  [2]  are  compared  with  chose  obtained  for  human  skin  in  vivo.  It  is 

seen  chat  Che  dielectric  behavior  of  human  skin  in  vivo  as  well  as 
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excised  skin  tissue  is  characterized  by  the  relaxation  of  free  water. 
Beside  the  expected  differences  between  excised  and  live  tissues 
regarding  water,  protein  and  salt  content,  it  should  be  noted  chat  the 
results  for  excised  skin  tissue  are  at  37#C  while  the  live  tissue 
results  are  at  the  ambient  temperature  (23°C).  The  difference  in 
temperature  could  accounc  for  differences  in  relaxation  time  between 
live  skin  and  excised  skin  tissue.  The  relaxation'  time  of  the  water  in 
live  skin  (9.205  ps)  is  somewhat  higher  than  the  relaxation  time  of  pure 
water  at  23®C  (8.73  ps).  This  may  be  due  to  hydrogen  bonding  between 
water  molecules  and  active  groups  of  other  types  of  molecules  [2]. 

Schwan  f  7  ]  has  previously  observed  that  the  low  frequency 
permittivity  (eg)  is  determined  by  the  total  solid  content  of  the 
tissue.  Assuming  dielectric  decrement  of  1  per  gm  solids  in  100  c.c. 
[7],  p  which  denotes  percent  solid  content  (gms/100  c.c.)  is  given  by 


where  Egw  is  the  low  frequency  permittivity  of  water.  From  [22]  for 

water  at  23°C  is  78.56.  So,  the  e  of  25.44  for  human  skin  in  vivo 

sw  -  ■  '  - 

corresponds  to  53.12  percent  solid  content  while  the  eg  of  42  for 
excised  skin  tissue  corresponds  to  36.56  percent  solid  content.  So,  the 
reduction  in  tg  for  human  skin  in  vivo  as  compared  with  excised  skin 
tissue  can  be  explained  on  the  basis  of  solid  content  of  the  skin 
tissue. 


8  L 


REFERENCES 


1.  0.  P.  Gandhi  and  A.  Riazi,  -Absorption  of  millimeter  waves  by  human 

beings  and  its  biological  implications,”  IEEE  Transactions  on 
Miarouave  Theory  and  Technique,  Vol.  34,  pp.  228-235,  1986. 

2.  H.  P.  Cook,  "The  dielectric  behavior  of  some  types  of  human  tissues 

at  microwave  frequencies,"  British  Journal  of  Applied  Physics,  Vol. 
2,  pp.  295-300,  1951. 

3.  T.  S.  England,  "Dielectric  properties  of  the  human  body  in  the 

microwave  region  of  the  spectrum,"  Nature,  Vol.  163,  pp.  487-488, 

1949. 

4.  T.  S.  England,  “Dielectric  properties  of  the  human  body  for  the 

wavelengths  in  the  1  -  10  cm  range,"  Nature,  Vol.  166,  pp.  480—481, 

1950. 

5.  H.  P.  Schwan  and  K.  Li,  "The  mechanism  of  absorption  of  ultrahigh 
frequency  electromagnetic  energy  in  tissues,  as  related  to  the 
problem  of  tolerance  dosage,"  IPE  Transactions  on  Medical 
Electronics,  Vol.  4,  pp.  45-49,  1956. 

6.  H.  P.  Schwan  and  K.  Li,  "Hazards  due  to  total  body  irradiation  by 
radar,"  Proceedings  of  IRE ,  Vol.  44,  pp.  1572-1581,  1956. 

7.H.  P.  Schwan,  "Electrical  properties  of  tissue  and  cell 
suspensions,"  Advances  in  Biological  and  Medical  Physics,  Vol.  5, 
pp.  147-209,  1957. 

8.  G.  L.  Hey-Shipton,  P.  A.  Matthews  and  J.  Mestay,  "The  complex 

permittivity  of  human  tissue  at  microwave  frequencies,"  Phys.,  Med. 
Biol.,  Vol.  27,  pp.  1067-1071,  1982. 

9.  D.  K.  Ghodgaonkar,  M.  F.  Iskander,  and  0.  P.  Gandhi,  "Complex 

permittivity  measurement  of  biological  tissues  and  human  skin  in 
vivo  at  millimeter  wavelengths,"  Bioelectramagnetics  Society 
Meeting,  June  16-20,  1985,  San  Francisco,  California. 

10.  J.  Galejs,  "Admittance  of  a  waveguide  radiating  into  stratified 

plasma,"  IEEE  Transactions  on  Antennas  and  Propagation,  Vol.  13, 

pp.  64-70,  1965. 

11.  C.  C.  Johnson  and  A.  W.  Guy,  "Nonionizing  electromagnetic  wave 

effects  in  biological  materials  and  systems,"  Proceedings  of  the 
IEEE,  Vol.  60,  pp.  692-718,  1972. 

12.  R.  Caldecott,  "The  generalized  multiprobe  ref lectoraeter  and  its 
application  to  automated  transmission  line  measurements,"  IEEE 
Transactions  on  Antennas  and  Propagation,  Vol.  21,  pp.  550-555, 
1973  . 


13.  D.  E.  duller,  "A  method  for  solving  algebraic  equations  using  an 
automated  computer,”  Matnemattca^  Tables  ana  Aias  to  computation , 
Vol.  10,  pp.  208-215,  1956. 

14.  N.  Marcuvitz,  Waveguide  Handbook ,  Me  Graw-Hill  Book  Company,  Inc., 
New  York,  1^51 . 

15.  M.  C.  Decreton  and  F.  E.  Cardiol,  "Simple  nondestructive  method  for 
the  measurement  of  complex  permittivity,"  IEEE  Transactions  on 
Instrumentation  and  Measurement ,  Vol.  23,  pp.  434-438,  1974. 

16.  W.  Culshah  and  M.  V.  Anderson,  "Measurement  of  permittivity  and 
dielectric  loss  with  a  millimeter  wave  Fabry-Parot  Interferometer," 
Proa.  Inst .  Ele.  Eng.,  Suppl.  23,  Vol.  109,  pp.  820-826,  1962. 

17.  J.  E.  Degenford,  "A  quasi-optic  technique  for  measuring  dielectric 
loss  tangent,  IEEE  Transactions  Instrumentation  and  Measurement, 
Vol.  17,  pp.  413-417,  1968. 

18.  S.  Roberts  and  A.  Von  Hippel,  "A  new  method  for  measuring 
dielectric  conscanc  and  loss  in  the  range  of  centimeter  waves,” 
Journal  of  Applied  Physics,  Vol.  17,  p.  610,  1946. 

19.  D.  K.  Ghodgaonkar,  0.  P.  Gandhi,  and  M.  F.  Iskander,  "Measurement 
of  complex  permittivities  of  biological  materials  in  the  frequency 
band  26.5  -  60  GHz,"  submitted  for  publication  to  IEEE  Transactions 
on  Microuave  Theory  and  Techniques. 

20.  M.  F.  Iskander  and  S.  S.  Stuchly,  "A  time-domain  technique  for 
measurement  of  the  dielectric  properties  of  biological  substances," 
IEEE  Transactions  on  Instrumentation  and  Measurement ,  Vol.  21,  pp. 
425-429,  1972. 

21.  R.  J.  Sheppard,  B.  P.  Jorden  and  E.  H.  Grant,  "Least-sauare 
analysis  of  complex  data  with  application  to  permittivity 
measurements,"  Journal  of  Physics,  D:  Applied  Physics,  Vol.  3,  pp. 
1759-1764,  1970. 

22.  E.  H.  Grant  and  R.  Shack,  "Complex  permittivity  measurements  at  8.6 
mm  wavelength  over  the  temperature  range  1  -  60*C,"  British  Journal 
of  Applied  Physics,  Vol.  18,  pp.  1807-1813,  1967. 


-  83  - 


Table  1. 


Debye  parameters  for  human  skin  In  vivo, 
and  excised  human  skin  tissue* 


Fig.  1.  Block  diagram  of  automated  measurement  system  for  the  frequency  band 
40  -  60  GHz.  Similar  measurement  system  has  been  developed  for 
26.5  -  40  GHz  band. 
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Fig.  3.  Flow  chart  for  in  vivo  dielectric  measurements. 


Fig.  7.  Emissivity  (percent)  for  the  human  skin 
in  vivo  (present  results)  and  the 
excised  skin  tissue. 


Complex  permittivity  values  for  human  skin 
In  vivo  (present  results)  and  excised 
skin  tissue. 
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